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ABSTRACT 


We study the clustering properties of metals in the intergalactic medium (IGM) 
as traced by 619 C IV and 81 Si IV absorption components with N > 10^^ cm“^ 
and 316 Mg II and 82 Fe II absorption components with N > 10^^’® cm“^ in 19 high 
signal-to-noise (60-100 per pixel), high resolution {R = 45000) quasar spectra. C iv 
and Si IV trace each other closely and their line-of-sight correlation functions ^{v) 
exhibit a steep decline at large separations and a flatter profile below « 150 km s“^, 
with a large overall bias. These features do not depend on absorber column densities, 
although there are hints that the overall amplitude of ^civ increases with time over 
the redshifts range detected (1.5-3). Carrying out a detailed smoothed particle hydro- 
dynamic simulation (2 x 320^, 57 Mpc^ comoving), we show that the C iv correlation 
function can not be reproduced by models in which the IGM metallicity is a constant 
or a local function of overdensity {Z oc A^/^). However, the properties of ^civ(ii) are 
generally consistent with a model in which metals are confined within bubbles with 
a typical radius Rg about sources of mass > Mg. We derive best-fit values of ~ 2 
comoving Mpc and Mg ss 10^^ M© at z = 3. Our lower redshift (0.5 — 2) measure¬ 
ments of the Mg II and Fe ii correlation functions also uncover a steep decline at large 
separations and a flatter profile at small separations, but the clustering is even higher 
than in the z =1.5-3 measurements, and the turn-over is shifted to somewhat smaller 
distances « 75 km s“^. Again these features do not change with column density, but 
there are hints that the amplitudes of ^Mgii('c) and ^Feii('c) increase with time. We 
describe an analytic “bubble” model for these species, which come from regions that 
are too compact to be accurately simulated numerically, deriving best-fit values of 
Rg w 2.4 Mpc and Mg « 10^^ Mq. Equally good analytic fits to all four species are 
found in a similarly biased high-redshift enrichment model in which metals are placed 
within 2.4 comoving Mpc of Mg ^ 3 x 10® sources at z = 7.5. 


1 INTRODUCTION 

Pollution is ubiquitous. Even in the tenuous intergalactic 
medium (IGM), quasar (QSO) absorption line studies have 
encountered heavy elements in all regions in which they were 
detectable (Tytler et al. 1995; Songaila & Cowie 1996). Such 
analyses were limited at first to somewhat overdense re¬ 
gions of space, traced by Lya clouds with column densities 
Ahi > ® cm“^. Here measurements of Nciv/Nm indi¬ 

cated that typically [C/H] ~ —2.5 at z ~ 3, with an order 
of magnitude scatter (Hellsten et al. 1997; Rauch, Haehnelt, 
& Steinmetz 1997). 

Pushing into more tenuous regions, statistical meth¬ 
ods have shown that unrecognized weak absorbers must be 


present in order to reproduce the global C iv optical depth 
(Ellison et al. 2000), and that a minimum IGM metallicity 
of approximately 3 x 10~^Zq was already in place at z = 5 
(Songaila 2001; hereafter SOI). While the filling factor of 
metals in such tenuous structures is an object of intense in¬ 
vestigation and debate (Schaye et al. 2000; Petitjean 2001; 
Bergeron et al. 2002; Carswell et al. 2002; Simcoe et al. 2002; 
Pettini et al. 2003; Schaye et al. 2003; Aracil et al. 2004) their 
very existence has profound cosmological implications. 

As the presence of metals increases the number of lines 
available for radiative cooling, even modest levels of enrich¬ 
ment can greatly enhance the cooling rate (e.g. Sutherland 
& Dopita 1993), which has the potential to accelerate the 
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formation of massive ( Xl, 10^^ M©) galaxies (e.g. Thacker, 
Scannapieco & Davis 2002). Furthermore, significant preen¬ 
richment is necessary to reproduce the abundances of G- 
dwarf stars in the Milky Way (e.g. van de Bergh 1962; 
Schmidt 1963) and nearby galaxies (e.g. Thomas, Greggio, 
& Bender 1999). 

Similarly, the violent events that propelled heavy ele¬ 
ments into the space between galaxies have important im¬ 
plications for the thermal and velocity structure of the IGM 
(e.g. Tegmark, Silk, & Evrard 1993; Gnedin & ©striker 1997; 
Gen & Bryan 2001). Outflows energetic enough to eject met¬ 
als from the potential wells of dwarf galaxies, for example, 
would have exerted strong feedback effects on nearby ob¬ 
jects (Thacker, Scannapieco, & Davis 2002). In this case the 
winds impinging on pre-virialized overdense regions would 
have been sufficiently powerful to strip the baryons from 
their associated dark matter, greatly reducing the num¬ 
ber of 10^° M© galaxies formed (Scannapieco, Ferrara, 
& Broadhurst 2000; Sigward, Ferrara, & Scannapieco 2005). 

Yet despite their many consequences, the details of how 
metals came to enrich the IGM are unclear. While numerous 
starburst-driven outflows have been observed at « = 3 (Pet- 
tini et al. 2001) and in lensed galaxies at 4 2 5 (Frye, 

Broadhurst, & Benitez 2002), it is unclear whether these ob¬ 
jects are responsible for the majority of cosmological enrich¬ 
ment. In fact a variety of theoretical arguments suggest that 
such galaxies represent only the tail end of a larger popula¬ 
tion of smaller “pre-galactic” starbursts that mostly formed 
at much higher redshifts (Madau, Ferrara, & Rees 2001; 
Scannapieco, Ferrara, & Madau 2002). On the other hand, 
active galactic nuclei are observed to host massive outflows 
(Begelman, Blandford, & Rees 1984; Weyman 1997), whose 
contribution from less luminous objects at intermediate red- 
shifts remains unknown (e.g. Fan et al. 2001). The impact 
of such lower-redshift events on the IGM is also hinted at by 
the “stirring” of C iv systems observed in studies of lensed 
QSO pairs (Rauch, Sargent, & Barlow 2001). Finally, a num¬ 
ber of theoretical studies suggest that primordial, metal- 
free stars may have been very massive (e.g. Bromm et al. 
2001; Schneider et al. 2002), resulting in a large number of 
tremendously powerful pair-production supernovae, which 
distributed metals into the IGM at extremely early redshifts 
15 (Bromm 2003; Norman, O’Shea, & Paschos 2004). 

While perhaps the main feature shared by such scenar¬ 
ios is their dependence on a poorly-understood population 
of presently undetectable objects, this assessment paints an 
overly bleak picture. Regardless of which objects enriched 
the IGM, it is clear that they must have formed in the dens¬ 
est regions of space, regions that are far more clustered than 
the overall dark matter distribution. Furthermore this “ge¬ 
ometrical biasing” is a systematic function of the masses of 
these structures, an effect that has been well studied analyt¬ 
ically and numerically (e.g. Kaiser 1984; Jing 1999). Thus 
the observed large-scale clustering of metal absorbers en¬ 
codes valuable information about the masses of the objects 
from which they were ejected. Likewise, as the maximal ex¬ 
tent of each enriched region is directly dependent on the 
velocity at which the metals were dispersed, measurements 
of the small-scale clustering of these absorbers are likely to 
constrain the energetics of their sources. 

Previous studies of the two-point correlation function of 
C IV components have shown that they cluster strongly on 


velocity scales up to 500 km s“^ (Sargent et al. 1980, Stei- 
del 1990, Petitjean & Bergeron 1994, Rauch et al. 1996). It 
has often been suggested that the clustering signal reflects a 
combination of (i) relative motions of clouds within a galac¬ 
tic halo and (ii) galaxy clustering. More recently Boksen- 
berg, Sargent, & Rauch (2003; hereafter BSR03) have gath¬ 
ered a sample of 908 C iv absorber components clumped 
into 199 systems in the redshift range 1.6 < 2 < 4.4 identi¬ 
fied in the Keck spectra of nine QSOs. They conclude that 
most of the signal is due to the clustering of components 
within each system, where a system is defined as a set of 
components that is “well-separated” from its neighbours as 
identified by the observer. In this case almost all the sys¬ 
tems extend less than 300 km s“^ and most extend less 
than 150 km s“^. They did not observe clustering between 
systems on the larger scales expected for galaxy clustering, 
although they concluded from their measurements of com¬ 
ponent clustering and ionisation balance that each system 
was closely associated with a galaxy. 

In Pichon et al. (2003; hereafter Paper I) we used 643 
C IV and 104 Si iv absorber components, measured by an 
automated procedure in 19 high signal-to-noise quasar spec¬ 
tra, to place strong constraints on the number and spatial 
distribution of intergalactic metals at intermediate redshifts 
(2 < 2 < 3). In this work, we showed that the correlation 
functions of intergalactic C iv and Si iv could be understood 
in terms of the clustering of metal bubbles of a typical co¬ 
moving radius Rs around sources whose biased clustering 
was parameterized by a mass Mg. A similar picture was also 
put forward in BSR03, but in our case significant large-scale 
clustering, similar to that seen in galaxies, was observed. 

In this paper we extend the analysis in Paper I in three 
important ways. First we carry out a more detailed study of 
the physical properties of C iv and Si iv absorbers and the 
relationship between local quantities and the overall spa¬ 
tial distribution. Second, we carry out a similar analysis 
of Mg II and Fe ii absorbers in our observational sample, 
which probe the IGM in a somewhat lower redshift range. 
Finally, we replace our dark-matter only modeling of Paper 
I with a full-scale smoothed-particle hydrodynamical sim¬ 
ulation. We then generate simulated metal-line spectra by 
painting bubbles of metals directly onto the gas distribu¬ 
tion at 2 > 2. By analyzing the resulting spectra with the 
same automated procedure applied to the measured data¬ 
set, we are able to place our models and observations on the 
same footing, drawing important constraints on the sources 
of metals. Motivated by measurements of the cosmic mi¬ 
crowave background, the number abundance of galaxy clus¬ 
ters, and high-redshift supernovae (e.g. Spergel et al. 2003; 
Eke et al. 1996; Perlmutter et al. 1999) we adopt cosmologi¬ 
cal parameters oih — 0.7, Q.m = 0.3, Da = 0.7, Db = 0.044, 
throughout this investigation where h is the Hubble Con¬ 
stant in units of 100 km s“^ Mpc“^ and Dm, Da, and fli, 
are the total matter, vacuum, and baryonic densities in units 
of the critical density, Pcrit. 

The structure of this work is as follows. In §2 we sum¬ 
marize the properties of our data set and reduction meth¬ 
ods. In §3 we present the number densities of C iv. Si iv. 
Mg II, and Fe ll, and estimate the cosmological densities of 
these species. In §4 we study the spatial clustering of these 
species and how it is related to local quantities such as col¬ 
umn density and abundance ratios. In §5 we describe our 
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numerical model for the distribution of neutral hydrogen in 
the IGM and compare it with observations. In §6 we ex¬ 
tend our method to include various models of cosmological 
enrichment and in §7 we compare these models to the ob¬ 
served distribution of C iv to derive constraints on the sizes 
and properties of sources of cosmological metals. In §8 we 
discuss an analytic model that is particularly suitable for 
comparisons with the distribution of Mg ll and Fe ll, as nu¬ 
merical analyses of these species are beyond the capabilities 
of our simulation. Conclusions are given in §9. 

2 DATA SET AND ANALYSIS METHODS 

2.1 Data and Reduction 

The ESO Large Programme “The Cosmic Evolution of the 
ICM” was devised to provide a homogeneous sample of QSO 
sight-lines suitable for studying the Lyman-a forest in the 
redshift range 1.7—4.5. High resolution (B « 45000), high 
signal-to-noise (60-100 per pixel) spectra were taken over 
the wavelength ranges 3100-5400 and 5450-10000 A, using 
the UVES spectrograph on the Very Large Telescope (VLT). 
Emphasis was given to lower redshifts to take advantage of 
the very good sensitivity of UVES in the blue and the fact 
that the Lyman-a forest is less blended. The distribution of 
redshifts, and the resulting coverage of various metal line 
absorbers are given in Table 1. In all cases we consider only 
metal absorption lines redward of the Lya forest, to avoid 
the extensive blending in this region, and blueward of 8110 
A, to avoid contamination from sky lines. The regions be¬ 
tween 5750-5830 A, 6275-6323 A, 6864-6968 A, 7165-7324 A, 
and 7591-7721 A were also excluded from our sample due 
to sky-line contamination. The C iv, Si iv, Mg ii, & Fe ii 
metal lines discussed in this paper were well-detected over 
the redshift ranges of 1.5-3.0, 1.8-3.0, 0.4-1.8, and 0.5-2.4 
respectively. 

Observations were performed in service mode over a 
period of two years. The data were reduced using the 
UVES context of the ESO MIDAS data reduction pack¬ 
age, applying the optimal extraction method, and follow¬ 
ing the pipeline reduction step by step. The extraction slit 
length was adjusted to optimize sky-background subtrac¬ 
tion. While this procedure systematically underestimates 
the sky-background signal, the final accuracy is better than 
1%. Wavelengths were corrected to vacuum-heliocentric val¬ 
ues and individual ID spectra were combined using a sliding 
window and weighting the signal by the total errors in each 
pixel. 

The underlying emission spectrum of each quasar was 
estimated using an automated iterative procedure that min¬ 
imizes the sum of a regularisation term and a term that 
was computed from the difference between the quasar spec¬ 
trum and the continuum estimated during the previous iter¬ 
ation. Finally the spectrum was divided by this continuum, 
leaving only the information relative to absorption features. 

2.2 Metal Line Identification 

Metal-line absorbers were identified using an automated 
two-step procedure. For each species that has multiple tran¬ 
sitions, we estimated the minimal flux compatible with the 


data for all pixels of the spectrum. This was done by first 
finding the pixels associated with the transition wavelengths 
Wi of a given species and then taking the maximum of the 
flux values in these pixels, scaled by Wifi, where fi is the 
oscillator strength associated with each of the transitions. 

A standard detection threshold was then applied to 
these spectra, such that only absorption features with equiv¬ 
alent widths (EWs) larger than 5 times the noise rms were 
accepted, giving a first list of possible identifications. This 
list was cleaned, using the similarity of the profiles of the 
transitions of a species and applying simple physical crite¬ 
ria that correlate the detection of two different species. For 
instance, one criterion implies that the detection of a Si iv 
system at a given redshift should be associated with the 
detection of a C iv system. 

Next, each system was fitted with Voigt profiles, taking 
care of their identification and possible blends with other 
systems. The first guess and the final Voigt profile decompo¬ 
sition were carried out using the VPFIT software (Carswell 
et al. 1987). Our decomposition of saturated systems is con¬ 
servative, in that it introduces additional unsaturated com¬ 
ponents only if there is some structure in the 1551 line that 
reveals their presence. This fitting procedure is described in 
detail in Aracil et al. (2005) and has been tested on sim¬ 
ulated spectra, doing well for all components with realistic 
values of N and b. 

Finally we applied a set of five cuts to the automated 
list generated by VPFIT: log N [cm“^] > 12 for C iv and 
Si IV and log N [cm~^] > 11.5 for Mg ll and Fe ll due to the 
detection limit of our procedure, 6 > 3 km s“^ to avoid false 
detections due to noise spikes, log A[cm“^] < 16 to remove 
very badly saturated components, and fe < 45 km s“^ to 
avoid false detections due to errors in continuum fitting. For 
the analyses presented here, we removed all associated com¬ 
ponents within 5000 km s“^ of the quasar redshifts. These 
cuts resulted in a final data set of 619 C iv (1548 A, 1551 
A), 81 Si IV (1394 A, 1403 A), 316 Mg ii (2796 A, 2803 A), 
and 82 Fe ll (2344 A, 2473 A, 2382 A) components, drawn 
from 688, 102, 320, and 88 components respectively, if we 
include the associates. These numbers differ slightly from 
those presented in Paper I due to further refinements in our 
detection procedure. 


3 NUMBER DENSITIES 

We first used our sample to compute the column den¬ 
sity distribution function, f(N), again working in the 
above assumed cosmology. Following Tytler (1987), /(A) 
is defined as the number of absorbing components per 
unit column density and per unit redshift path, dX. 
In this paper, we adopt a definition of dX = (l-b 
2 )^ [Da-I-D m(l-I- 2 )^] dz such that at all redshifts 
/(X) does not evolve for a population whose physical size 
and comoving space density are constant. Note that this 
definition is slightly different from that used in Paper I 
and in SOI, namely dX' = (1 -|- zf^^dz, although when 
2 > (DA/flm)^/® - 1 = 0.32, as is appropriate for our sam¬ 
ple, dX' can be very closely approximated as for 

comparison with previous analyses. 

In FigureQwe plot /(X) for both C iv and Si iv com¬ 
ponents, as was presented in Paper I. The mean redshifts 
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Name 

^em 

Eorest 

C IV 

Coverage 

Si IV 

Mg II 

Ee II 

PKS 2126-158 

3.280 

2.61-3.28 

2.36-3.28 

2.74-3.28 

0.85-1.85 

1.03-2.42 

Q 0420-388 

3.117 

2.47-3.12 

2.23-3.12 

2.59-3.12 

0.79-1.85 

0.95-2.42 

HE 0940-1050 

3.084 

2.45-3.08 

2.21-3.08 

2.56-3.08 

0.77-1.85 

0.93-2.42 

HE 2347-4342 

2.871 

2.27-2.87 

2.04-2.87 

2.38-2.87 

0.68-1.85 

0.83-2.42 

HE 0151-4326 

2.789 

2.20-2.79 

1.97-2.79 

2.31-2.79 

0.64-1.85 

0.79-2.42 

Q 0002-422 

2.767 

2.18-2.77 

1.96-2.77 

2.29-2.77 

0.64-1.85 

0.78-2.42 

PKS 0329-255 

2.703 

2.13-2.70 

1.91-2.70 

2.23-2.70 

0.61-1.85 

0.75-2.42 

Q 0453-423 

2.658 

2.09-2.66 

1.87-2.66 

2.19-2.66 

0.59-1.85 

0.73-2.42 

HE 1347-2457 

2.611 

2.05-2.61 

1.83-2.61 

2.15-2.61 

0.57-1.85 

0.70-2.42 

HE 1158-184.3 

2.449 

1.91-2.45 

1.71-2.45 

2.01-2.45 

0.50-1.85 

0.63-2.42 

Q 0329-385 

2.435 

1.90-2.44 

1.70-2.44 

2.00-2.44 

0.49-1.85 

0.62-2.42 

HE 2217-2818 

2.414 

1.88-2.41 

1.68-2.41 

1.98-2.41 

0.48-1.85 

0.61-2.41 

Q 1122-1328 

2.410 

1.87-2.41 

1.68-2.41 

1.98-2.41 

0.39-1.85 

0.61-2.41 

Q 0109-3518 

2.404 

1.87-2.40 

1.67-2.40 

1.97-2.40 

0.48-1.85 

0.61-2.40 

HE 0001-2340 

2.263 

1.75-2.26 

1.56-2.26 

1.84-2.26 

0.42-1.85 

0.54-2.26 

PKS 0237-23 

2.222 

1.72-2.22 

1.53-2.22 

1.81-2.22 

0.40-1.85 

0.53-2.22 

PKS 1448-232 

2.220 

1.72-2.22 

1.53-2.22 

1.81-2.22 

0.40-1.85 

0.52-2.22 

Q 0122-380 

2.190 

1.70-2.19 

1.50-2.19 

1.78-2.19 

0.38-1.85 

0.51-2.19 

HE 1341-1020 

2.135 

1.65-2.14 

1.46-2.14 

1.74-2.14 

0.36-1.85 

0.48-2.14 


Table 1. List of lines of sight. Here Zem is the quasar redshift, while Lya forest is used only redward of the Ly/3 transition at 1025.7 A, 
and metal absorption lines are used only redward of the Lyo forest and blueward of 8130 A. 
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Figure 1. Column density distribution of C iv (upper panel) 
and Si iv (lower panel) absorption components. In each panel, 
components are divided into two redshift bins: 1.5 < ^ < 2.3 
(squares) and 2.3 < 2 < 3.1 (triangles). Column density bins are 
IqO.Sjy cm~^ wide and error bars in this and all further plots are 
l-(7. The dashed line is the power law fit measured in SOI. Finally 
the small crosses are the full set of C iv and Si iv components 
identified by BSR03, with our cuts imposed. 


of C IV and Si iv in our sample were 2.16 and 2.38, re¬ 
spectively, and so in this plot we divide the data into two 
redshift bins from 1.5 < 2 < 2.3 and 2.3 < 2 < 3.1. Both 
species are consistent with a lack of redshift evolution, as 
found by previous lower-resolution studies of C iv and Si iv 
(SOI; Pettini et al. 2003), and pixel-by-pixel analyses of in- 
tergalactic C iv (Schaye et al. 2003). The overall density 
distribution of C iv is also consistent with a power-law of 
the form f{N) = BN~°‘ with a = 1.8 and logjQ / = —12.7 
at 10^® cm~^ as fitted by SOI. Finally, we compare our re¬ 
sults with the dataset collected in BSR03 from nine QSO 
spectra with a S/N ~ 50 per pixel. Here and below we use 
the full data set taken by BSR03, to which we apply ex¬ 
actly the same cuts as we do to our data. For components 
with columns « 10^^ cm“^ these data sets are quite sim¬ 
ilar. However, a significant difference between this sample 
and our own is the fit to the saturated C iv components 
with log(A^civ) > 14. These have been decomposed into a 
large number of smaller log(A^civ) < 12.5 systems in the 
BSR03 analysis, while our decomposition only introduces 
additional unsaturated components if there is structure in 
the 1551 line. Extrapolating the results of Songaila (2001) 
to column depths below 10^® cm“^ also yields a distribution 
similar to ours. 

While fewer in total, the Si iv components in the lower 
panel are also consistent with a lack of evolution, following a 
similar power law with a lower overall magnitude. Note that 
in this figure the error bars are purely statistical, estimated 
as one over the square-root of the number of components 
in each bin. Again, for comparison, we include the number 
densities computed from the full BSR03 sample, with our 
cuts applied. While this comparison is noisier, the overall 
trends are the same: at 10^® cm“^ the number densities are 
similar, while saturated components are decomposed into a 
larger number of smaller systems in the BSR03 data set. 











^ g I I I I I I I I I I I I I I I M I I I I I 

11 12 13 14 15 


log N/cm ^ 


Figure 2. Column density distributions of Mg ii (upper panel) 
and Fe ii (lower panel) absorption components. In each panel, 
components are divided into two redshift bins: 0.4 < 2 < 1.15 
(squares) and 1.15 < 2 < 1.9 (triangles). As in Fig. 1, the column 
density bins are 10’^ ®A'' cm“^ wide. The dashed lines correspond 
to the power-law fits described in the text, and in the lower panel 
we also include / values when the large 2 = 0.836 system in 
Q0002-422 is removed (see text). 


In Figure[2lwe plot f{N) for both Mg ll and Fe ll, now 
going down to a minimum column density of cm“^ 

which corresponds to roughly the same optical depth as 10^^ 
cm“^ for C IV and Si iv. For Mg ll and Fe ll the relevant 
doublets are at substantially longer restframe wavelengths, 
and therefore our UVES detections primarily occur at lower 
redshifts. Thus the mean redshifts of Mg ll and Fe ll are 
only 1.05, and 1.38, and we divide our data into bins from 
0.4 < z < 1.15 and 1.15 < z < 1.9. These lines arise in 
lower ionization gas and are often thought of as tracers of 
quiescent clouds, probably associated with galaxies (e.g. Pe- 
titjean & Bergeron 1990; Churchill et al. 1999; Churchill, 
Vogt, & Charlton 2003). 

Like its higher-ionization counterparts. Mg ll is consis¬ 
tent with a lack of evolution in number densities over the 
observed redshift range. In the Fe ii case, however, a sig¬ 
nificant excess of intermediate column density components 
is found at lower redshifts. A closer inspection of the data 
indicates that this feature is caused by a single large system 
in Q0002-422, at 2 = 0.836, which spans over 560 km s“^. 
The removal of this system results in the third set of points 
in the lower panel of Fig. [3 which are consistent with the 
higher-redshift values. The large impact of this system in 
our measurements suggests that simple y/N estimates may 
somewhat underpredict the statistical error on our measure¬ 
ment. This hints at strong clustering between Fe ll compo¬ 
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nents, which is in fact measured, as we discuss in detail 
below. 

Statistical fluctuations aside, the overall density distri¬ 
butions of Mg II and Fe ii are largely consistent with the 
power law fits obtained from previous measurements, apart 
from showing only a weak deviation in the lowest Npe u bin, 
probably due to incompleteness. In this case, the dashed- 
lined fits in Fig 121 are f(JV) = BN~°‘ with a = 1.6 and 
logio / = —13.2 at 10^^ cm“^ for Mg ll and a = 1.7 and 
logio / = —13.4 at 10^® cm“^ for Fe ii. While some flatten¬ 
ing of /Mgii(V) at even higher columns is necessary to match 
observations at column densities > 10^®’® cm~^ (Prochter, 
Prochaska, & Buries 2004), for the column densities in our 
sample our measured slopes are identical with those deter¬ 
mined by previous studies. In particular our q fits match 
those of Churchill, Vogt, & Charlton (2003), although our 
B values are different as these authors did not attempt to 
normalize their results by the total redshift path observed. 

In summary, our automatic identification procedure 
produces a set of components whose column density dis¬ 
tributions are consistent with previous measurements, com¬ 
plete to N 10^^ cm“^ for C iv and Si iv, and complete 
to N 10^^'® cm“^ for Mg ii and Fe ii. No evolution in / 
is seen for any species over the full redshift range probed, 
indicating that the majority of IGM enrichment is likely to 
have occurred before the redshifts observed in our sample. 

Finally, our number densities allow us to compute the 
total cosmological densities of each of the detected species. 
Following SOI, we express these in terms of a mass fraction 
relative to the critical density, which can be computed as 


fli< 


Horriiou ^ ^ ^ 

cpcrit AX AX 


( 1 ) 


where Hq is the Hubble constant, mion is the mass of the 
given ion, A is its atomic number, and AX is the total red¬ 
shift path over which it is measured. The results of this anal¬ 
ysis are given in Table 2. Note that these values are species 
densities, and no ionization corrections have been applied to 
estimate the corresponding element densities. Again these 
values are broadly consistent with previous measurements, 
although there is a significant scatter due to the fact that 
most of the material lies in the largest, rarest components. 
Thus previous studies have found flciv values as disparate 
as 6.8 X 10~® at 2 = 2.5 (SOI), 3.8 ± 0.7 x 10~® (BRS03), 
and between 3.5 x 10“® and 7.9 x 10“® depending on the 
method of analysis (Simcoe, Sargent, & Rauch 2004). 


4 SPATIAL DISTRIBUTION 
4.1 Carbon IV and Silicon IV 

Having constructed a sample of well-identified metal absorp¬ 
tion components, we then computed their two-point correla¬ 
tion function in redshift space, ^(v). This quantity was pre¬ 
viously studied in Rauch et al. (1996) who noted a marked 
similarity between ^(v) of C iv and Mg ii, in BSR03, who 
carried out a two-Gaussian fit (see also Petitjean & Bergeron 
1990, 1994), and in Paper I. For each quasar, we computed 
a histogram of all velocity separations and divided by the 
number expected for a random distribution. Formally, the 
correlation function for a QSO £ is 
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Species 

(^> 

log(N/cm-2) 

n 

An 

C IV 

2.2 

12-16 

7.54 X 10“® 

±2.16 X 10“® 

Si IV 

2.4 

12-16 

6.00 X 10-® 

±1.21 X 10“® 

Mg II 

1.1 

11.5-16 

5.95 X 10-® 

±2.23 X 10“® 

Fe II 

1.4 

11.5-16 

1.87 X 10“® 

±0.36 X 10-® 


Table 2. Cosmological densities of detected species. 


+ 1 



( 2 ) 


where is the number of pairs separated by a velocity 
difference corresponding to a bin fe, and (wfc) is the average 
number of such pairs that would be found in the redshift 
interval covered by QSO I, given a random distribution of 
redshifts with an overall density equal to the mean density 
in the sample. Alternatively, we may consider all QSOs at 
once and compute 


+ 1 — 




(3) 


or equivalently 

+ 1 = + 1] with wi= ^ (4) 


that is weighting the correlation found for each QSO by 
the number of random pairs that are expected given the 
redshift coverage of that QSO. The statistical variance in 
this measurement is given by 

(5) 

e 


where is the variance associated with bin k of quasar 
i. In Paper I, we estimated this quantity according to the 
usual formula 


a 


2,1 

k 



( 6 ) 


which gives the Poisson error in our measurement. In the re¬ 
sults presented here, however, we adopt a more conservative 
approach, and also include the additional scatter caused by 
the finite sample size used to construct the correlation func¬ 
tion (Mo, Jing, & Borner 1992). In this case 


(7 


2 ,e 

k 



, A 

rik + 4: 


Ml!' 

m ’ 


(7) 


where is the total number of components detected in 
QSO £. Note that the presence of this additional scatter 
highlights the strength of our high signal-to-noise data set, 
as it allows us to work in the limit in which the number of 
C IV components detected in each quasar is large. 

The resulting correlation functions are shown in Fig. 
m again split into two redshift bins. Interestingly, in the 
better-measured C iv case, there are hints that the z < 
2.3 correlation function may be enhanced with respect to 
the high-redshift one. Furthermore, this growth is consistent 
with a population of absorbers that “passively” evolves by 



100 1000 
velocity (km/s) 


Figure 3. Two point correlation function of C iv (upper panel) 
and Si iv (lower panel) absorption components. In each panel 
the components have been divided into two redshifts bins, with 
symbols as in Fig. 1. The upper panel also includes a number of 
comparisons with previous measurements. In particular: the lower 
set of crosses corresponds to the full set of components defined 
in BSR03, normalizing each QSO individually; the upper set of 
crosses corresponds to imposing a column density cut of 10*^^ 
cm“^, normalizing each QSO by its expected number of pairs; 
and the solid lines correspond to dividing the BSR03 data into a 
subsamples with z >Z.l (lower line) and 2 : < 3.1 (upper line), as 
described in the text. In the lower panel, the crosses corresponds 
to imposing a column density cut of 10^^ cm“^ on the Si iv 
components observed in BSR03 and normalizing each QSO by its 
expected number of pairs. 


following the motion of the IGM during the formation of 
structure, as we discuss in further detail in §8. 

In the upper panel of this figure we also plot correla¬ 
tion functions computed from the sample defined in BSR03, 
which is drawn from the spectra of nine QSOs with a mean 
redshift of 3.1 and a signal-to-noise per pixel « 50. In 
this case we show results obtained both from using the full 
data set, normalizing each quasar individually (as was car¬ 
ried out in BSR03), and from imposing a lower cut-off at 
Aciv.min = 10^^ cm“^, normalizing each quasar by the ex- 
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100 1000 


velocity (km/s) 

Figure 4. Dependence of the C iv correlation function on column 
density threshold. Upper panel: Effect of applying a cut on the 
maximum column density of C iv components used to calculate 
?Civ(r’)- In all cases Vciv.min = 10^^ cm“^. Lower panel: Effect 
of applying a cut on the minimum C iv column density, with 
IVciv.max fixed at 10^® cm“^. 


pected number of pairs (as was carried out in our analysis). 
In both cases the resulting iv values are similar and some¬ 
what lower in amplitude than our measurements. Rauch et 
al. (1996) similarly have found a lower amplitude. Dividing 
the BSR03 data into a z < 3.1 bin with a mean redshift of 
2.5 and a. z > 3.1 bin with a mean redshift of 3.6 resulted 
in correlation functions given by the solid curves (again cal¬ 
culated according to our method). Furthermore the ampli¬ 
tude of the a = 2.5 BSR03 correlation function is similar 
to our measurements, which are drawn from a sample with 
a mean redshift of 2.3. However, the higher-redshift curve 
is substantially lower, again indicating that iv is likely to 
evolve with redshift. This was also suggested by the analysis 
in BSR03 Figure 14, although they point out the changing 
ionizing background may also be an issue. Finally we note 
that while the BSR03 sample shows a relative lack of com¬ 
ponents at « 500 km/s. This is very near the C iv doublet 
separation. 

Moving to the bottom panel, we see that the overall 
shape and amplitude of the C iv and Si iv correlation func¬ 
tions are similar and are consistent to within the Si iv mea¬ 
surement errors, as was discussed in Paper I. Both functions 
exhibit a steep decline at large separations and a flatter pro¬ 
file at small separations, with an elbow occurring at « 150 
km s”^. Both functions are also consistent with the corre¬ 
lation one obtains from the full BSR03 Si iv sample, after 
applying our cuts. Finally, as was noted in Paper I, there 
is a weak low redshift feature at « 500 km s“^ in ^civ, the 
origin of which we explore in §4.2. 



Figure 5. The average log(Afc iv/Vsiiv) ratio for C iv compo¬ 
nents contributing to the correlation function at various separa¬ 
tions. At each separation the dashed error bars are the statistical 
errors, while the solid error bars are the intrinsic scatter. 

In Fig. 131 we study the dependence of the C iv spatial 
distribution on column density, by computing the correla¬ 
tion function over the full redshift range but selecting com¬ 
ponents within a fixed range of column density. In the upper 
panel of this figure we apply a cut on the maximum column 
density component, while holding the minimum Nciv fixed 
at our detection limit of 10^^ cm“^. Apart from a weak 
shift in the 500-630 km s“^ bin, ^civ(«) remains practically 
unchanged by this threshold. As the majority of the de¬ 
tected components are relatively weak, this indicates that 
our signal is determined by the bulk of the components in 
our sample, rather than by properties of individual strong 
absorbers. 

The results of a more drastic test are shown in the 
lower panel of this figure. Here we hold Nc iv,max fixed at 
10^® cm~^ and apply a cut on the minimum column den¬ 
sity, which greatly reduces the number of components in the 
sample. Nevertheless moving from Nciv,min = 10^^ cm“^ 
to Nciv.min = cm“^ results only in a very weak en¬ 

hancement of ^civiv) at small separations, while the rest 
of the correlation function remains unchanged. Thus, unlike 
Lya absorption systems (Cristiani et al. 1997), the corre¬ 
lation of C IV does not depend strongly on absorption col¬ 
umn densities. Instead, the spatial distribution seems to be 
a global property of the population of C iv components. 

A question that immediately arises is whether the fea¬ 
tures observed in the C iv and Si iv correlation functions are 
intrinsic to the underlying distribution of metals, or perhaps 
arise from variations in the UV background at somewhat 
shorter wavelengths. In fact, analyses of the He ii distribu¬ 
tion due to ionization by 54.4 eV photons suggest that the 
second reionization of hydrogen may still have been quite 
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patchy at « = 2.3 — 2.9 (Shull et al. 2004), with He ll found 
preferentially in “void” regions where HI is weak or unde¬ 
tected. 

On the other hand, the ionization potentials of C ill and 
Si III are 47.5 eV and 33.5 eV, respectively, somewhat lower 
than that of He ii, but well beyond the ionization poten¬ 
tial of hydrogen. Thus if the suggested patchiness of He ll is 
due primarily to changes in the IGM opacity at wavelengths 
shortward of 54.4 eV, then the distribution of C iv and Si iv 
is likely to more closely trace the underlying distribution of 
metals. If He ll inhomogeneities exist and are caused by a 
sparsity of hard sources, however, it is possible that back¬ 
ground variations may also play a role in the distribution of 
triply ionized regions of carbon and silicon. 

As the ionization potentials of C iii and Si iii differ by 
12 eV, each is sensitive to a slightly different range of UV 
photons. Thus if the features seen in FigureOwere produced 
by changes in the ionizing background, one might expect 
to see systematic changes in the ratio of these species as a 
function of separation. As a simple test of this possibility, 
we considered the average \og{Nciv/Nsiiv) as a function of 
separation. In order to make the sample included in this 
average as large as possible we computed this as 


log( 


Nci- 

Nsii 


) = 


Eijgbin, log (fel) ^(5 - - ^i\) 




1 X 0(5 — \ve — Wi|) 


( 8 ) 


where 6{v) is the Heaveside step function, i and j are indexes 
of C IV components, I is an index over all Si iv components, 
and fc is a given bin in velocity separation used to calcu¬ 
late the correlation function. In other words, for each bin in 
the correlation function, we average log(A^‘^ over 

all C IV components i that are found at the appropriate 
separation from another C iv component j and within 5 km 
s“^ of a Si IV component 1. 

The results of this analysis are found in Figure 
which shows no correlation between separation and 
species abundances. Furthermore, our average value of 
log(A'^"'/A'®"’') « 0.7 is similar to that seen in previous 
analyses of 10^^ cm~^ < Nciv < 3 x 10^“^ cm“^ absorbers 
(Kim, Cristiani, & D’odorico 2002; BSR03), as well as the 
weaker C iv and Si iv lines detected by Aguirre et al. (2003) 
using the pixel optical depth method. Thus there is nothing 
particularly unusual about the subset of absorbers selected 
by our procedure. Although this is clearly not an exhaustive 
test, it nevertheless suggests that the features in the correla¬ 
tion functions are not imprinted in a straightforward way by 
the UV background itself, and are more likely to be caused 
by the spatial distribution of metals. However, a much more 
detailed analysis is necessary to settle this issue definitively. 


4.2 Peculiar Systems at Low-Redshift 

The C IV correlation functions in Figs.|31andE]hint at a sec¬ 
ondary bump at large separations. It is important to try and 
understand if this comes from the presence of few peculiar 
systems or if this is a generic feature of the C iv distri¬ 
bution. To this end we computed the correlation function 
for different samples, each time excluding one of the lines of 



Figure 6. Top: Impact of peculiar systems on the C iv corre¬ 
lation function. Here the square points are computed from the 
full sample, the circles are computed excluding the sightline to¬ 
wards PKS 0237—23, and the stars are computed excluding both 
PKS 0237—23 and the two sub-DLA systems, as described in 
the text. Bottom: Comparison between $civ for the full sample, 
including associated absorbers (crosses), and excluding C iv com¬ 
ponents with 5000 km s“^ of the quasar redshifts (squares). 


sight, and discovered that the signal comes from three QSOs, 
namely, PKS 0237-23, HE 0001-2340 and Q 0122-380. 

The first of these has been long known to be very 
peculiar. Indeed, a huge C iv complex is seen toward 
PKS 0237—23 at eleven different redshifts over the range 
1.596—1.676 (more than 10,000 km s“^) with three main 
subcomplexes at Zabs = 1.596, 1.657 and 1.674 (Boroson et 
al. 1978, Sargent et al. 1988). Furthermore Folz et al. (1993) 
searched the field around PKS 0237—23 for other QSOs to 
provide background sources against which the presence of 
absorption at the same redshifts could be investigated. They 
concluded that the complex can be interpreted as a real spa¬ 
tial overdensity of absorbing clouds with a transverse size 
comparable to its extent along the line of sight, that is of 
the order of 30 Mpc. The correlation function without this 
line of sight is shown in Fig. El 

Two other lines of sight display peculiar systems. At 
Zahs = 2.1851 toward HE 0001 — 2340 there is a sub-DLA 
system and the associated C iv system is spread over 
Ri450 km s“^. It is therefore difficult to know if the struc¬ 
ture there is due to large scales or more probably to the 
internal structure of the halo associated with this high den¬ 
sity peak. At Zabs = 1.9743 toward Q 0122—380, there is a 
double strong system spread over more than 500 km s“^. It 
is again difficult to know whether these absorptions reflect 
internal motions of highly disturbed gas. 

After these are removed, the most significant excess at 
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Figure 7. Upper panel: Two point correlation function of Mg ii, 
divided into two redshift bins as in Figure 2. Lower panel: Two 
point correlation function of Fe ii, divided into the same redshift 
bins. 

large separations is found in the 500-630 km s“^ bin. This ve¬ 
locity difference corresponds to the difference in wavelengths 
of the C IV doublet itself. In fact, it is interesting to note 
that this bin is the only one that is significantly reduced by 
applying a cut to eliminate the larger A^civ components, as 
was seen in Fig. 0 

As a further test of large-separation correlations, we 
have also computed 5civ(u) including the associated sys¬ 
tems, found within 5000 km s“^ of the redshifts of the QSOs 
in this sample. This is compared with the C iv correlation 
function for our standard sample in the lower panel of Fig.|^ 
At all separations, ^civ remains unchanged, thus indicating 
that associated systems are not distributed in a particularly 
unusual way, and do not contribute any significant features 
to 5civ(n) at « 500 km s“^, or any other separation. 

4.3 Iron II and Magnesium II 

We now turn our attention to the distribution of lower- 
redshift metals, as traced by Mg ii and Fe ii. Splitting the 
data into two redshift ranges yields the line-of-sight correla¬ 
tion functions shown in Fig.0 where again we have included 
both the Poisson and sample-size errors in our estimate of 
the variances. Like their high-redshift counterparts. Mg ll 
and Fe ll are found to trace each other closely. Their corre¬ 
lation functions are both relatively shallow at small separa¬ 
tions and fall off more steeply at large separations. Also like 
^Civ(n), both 5Mgii(u) and ^Feii(u) exhibit slight enhance¬ 
ments at lower redshifts, although again these excesses fall 
within the errors. 

Next we examine the dependence of the Mg ll spatial 



100 1000 


velocity (km/s) 

Figure 8. Dependence of Mg ii correlation function on column 
density threshold. Upper panel: Effect of applying a cut on the 
maximum column density of Mg ii subcomponents used to calcu¬ 
late CMgii(r’). In all cases NMgii,min = 10^^ ® cm“^. Lower panel: 
Effect of applying a cut on the minimum Mg ii column density, 
with NMgii.max fixed at 10^® cm“^. 


distribution on column density. Removing the strongest ab¬ 
sorbers in our sample before calculating ^Mgiilu) results in 
the values plotted in the upper panel of Figure]^ As in the 
C IV case, the Mg ii correlation function is not dominated 
by the clustering of large components, but rather remains 
almost unchanged as a function of Nmax, even when it is re¬ 
duced to 10^^ ® cm“^, excluding over a third of the systems. 
Similarly, raising the minimum column density from 10^^ ® 
cm~^ to 10^^ ® cm“^ does not boost 5Mgii(u), even though 
this excludes « 2/3 of the sample. 

In FigureElwe compare the correlation functions of C iv 
and Si iv with those of Mg ll and Fe ll. Note that the mean 
redshift of these lower ionization species is « 1.2, while for 
C IV and Si iv Zmean ~ 2.3. Thus our sample contains very 
few objects in which all four species can be directly com¬ 
pared. Nevertheless, a comparison of their redshift-space 
correlation functions reveals a number of important paral¬ 
lels. While ^(v) of all species decline steeply at large sep¬ 
arations and exhibit a turn-over at smaller velocity differ¬ 
ences, the transition between these two regimes is pushed to 
slightly smaller separations in the Mg ll and Fe ll case, and 
the fall-off at higher densities is more abrupt. 

Interestingly, the features seen in this distribution can 
be inferred from the original fitting to the distribution of 
velocity separations of Mg ii absorbers by Petitjean & Berg¬ 
eron (1990), using a remarkably small number of systems. 
Their data were fit with the sum of two Gaussian distribu¬ 
tions with similar overall weights and velocities dispersions 
of (T„ = 80 km s“^ and (t„ = 390 km s“^, which the authors 
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Figure 9. Measured correlation function of all metal-line com¬ 
ponents. Points are measurements from our sample, while the 
solid line is the Mg ii fit by Churchill, Vogt, & Charlton (2003), 
arbitrarily normalized. The dashed line is the C iv correlation 
function, shifted upwards by a factor of 2.1 to provide a simple 
estimate of the impact of structure formation from 2 = 2.2 to 1.1 
on a fixed population of absorbers. 


interpreted as due to the kinematics of clouds bound within 
a given galaxy halo, and the kinematics of galaxies pairs, re¬ 
spectively. Working at higher spectral resolution and signal- 
to-noise, Churchill, Vogt, & Charlton (2003) also obtained 
a good two-component Gaussian fit to the two-point clus¬ 
tering function of Mg ll components, although they did not 
attempt to normalize this function to obtain ^Mgii(u) +1. In 
this case the best-fit values were = 54 km s“^ and (t„ = 
166 km s“^, where the relative amplitude of the narrow com¬ 
ponent was twice that of the broad component. This fit has 
been added to Figure adopting an arbitrary normaliza¬ 
tion. Although our dataset has an overall signal-to-noise that 
is higher than that of Churchill, Vogt, & Charlton (2003), 
and thus is more complete at lower column densities, their 
two-Gaussian model also provide a good match to our data 
at Av ^ 400 km/s. However, it falls short of the observed 
correlation at larger separations. 

To contrast the correlation functions in more detail, we 
have added a simple estimate of “passive” evolution to Figure 
El that is the evolution if the metals detected at 2 « 2.3 
as C IV absorbers were to move along with the formation 
of structure before appearing as Mg ll absorbers at 2 = 
1.2. To first approximation, the overall bias of such a metal 
tracer field would remained fixed, but its correlation function 
would be enhanced by a factor of Il^(1.2)/D^(2.3) = 2.1, 
where D{z) is the linear growth factor. Surprisingly, simply 
shifting 5c iv by a factor of 2.3 provides us with an accurate 
match for the Mg ll correlation function over a large range 
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velocity (km/s) 

Figure 10. Upper panel: Impact of linking together C iv compo¬ 
nents into systems. Lower panel: Impact of linking together Mg ii 
components into systems. Details described in the text. 

of separations, although it underpredicts the clustering of 
Mg II and Fe ii at smaller distances. This is dicussed in 
futher detail in §8. 

To facilitate future comparisons, in Table 3 we give the 
correlation function and errors for each of the four species 
averaged over our full sample. Note that the small number 
of Si IV and Fe ll components forces us to use a smaller 
number of bins to beat down the statistical noise in our 
measurements. 

Finally, we carry out a test to determine if the spa¬ 
tial distribution of metals as traced by 5civ(n) may be af¬ 
fected by our VPFIT decomposition into components. Pre¬ 
vious studies have attempted to trace the distribution of 
intergalactic metals by grouping together components into 
“systems,” which are likely to have a common physical ori¬ 
gin, and computing the correlation function of these systems 
(e.g. Petitjean & Bergeron 1990; BSR03). While typically, 
system identifications have been carried out by eye, here 
we attempt a more objective approach, which parallels the 
friends-of-friends technique (Davis et al. 1985) widely used 
for group-finding in cosmological simulations. In this case, 
we define a velocity linking-length (wunk) and group together 
all components whose separation from their nearest neighbor 
is less than nunk into a system at a redshift equal to the av¬ 
erage over all its components. Note that this procedure does 
not involve simply linking together pairs within uunk, but 
rather forms collections of many components, each within a 
linking length of its neighbors and grouped together into a 
single entity. It is therefore equivalent to partitioning a set of 
components into two systems whenever they are separated 
by a gap wider than uunk. 
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bin (km s 

?Civ 

CMg ii 

bin (km s 

^Si iv 

?Feii 

20-25 

41 ± 8 

170 ± 50 

20-30 

94 ± 52 

310 ± 150 

25-32 

66 ± 13 

170 ± 40 

30-43 

71 ± 30 

280 ± 130 

32-40 

59 ± 11 

240 ± 60 

4.3-65 

71 ± 36 

220 ± 100 

40-50 

40 ± 6 

140 ± 40 

65-100 

30 ± 17 

160 ± 70 

50-63 

49 ± 9 

145 ± 30 

100-140 

27 ± 13 

74 ± 34 

63-79 

35 ± 5 

155 ± 40 

140-200 

11 ± 6 

45 ± 19 

79-100 

30 ± 5 

96 ± 23 

200-300 

6.2 ± 4.2 

25 ± 11 

100-125 

26 ± 4 

98 ± 21 

300-450 

6.6 ± 3.5 

7.8 ± 4.1 

125-160 

20 ± 3 

64 ± 15 

450-670 

3.9 ± 2.6 

0.88 ± 0.58 

160-200 

14 ± 2 

42 ± 10 

670-1000 

0.8 ± 0.7 

Oil 

200-250 

7.3 ± 1.2 

35 ± 8 




250-320 

5.4 ± 1.0 

20 ± 5 




329-400 

3.6 ± 0.7 

12 ± 3 




400-500 

3.5 ± 0.7 

6.9 ± 2.3 




500-630 

3.6 ± 0.8 

1.7 ± 0.8 




630-790 

1.6 ± 0.4 

4.2 ± 2.4 




790-1000 

0.98 ± 0.32 

2.2 ± 1.3 





Table 3. Summary of measured metal-component correlation functions. Note that there is likely to be significant redshift evolution of 
these functions. The mean redshifts of C iv and Si iv are ~ 2.3. The mean redshifts of Mg ii and Fe ii are ~ 1.15. 


In the upper panel of Figure HTil we plot ^civ(n) com¬ 
puted for the resulting C iv systems, for three different 
choices of funk. In all cases, within our measurement er¬ 
rors, combining components into systems has no appreciable 
impact at separations much larger than the linking length. 
Thus while BSR03 report a lack of clustering of systems as 
identified by eye, we are unable to reproduce this behavior 
with our automatic method. Perhaps this is not surprising, 
as the clustering of ^civ(u) is very strong, and thus many 
pairs of “systems” are likely to be closely spaced and easily 
tagged as a single object. However, our results show that 
fixing a pre-specified definition of systems does not remove 
large-scale correlations in this way. 

In the lower panel of Figure uni we see that grouping 
Mg II components into systems has no clear impact at larger 
separations if uunk = 25, 50, and while there are hints of 
weak larger-scale damping if uunk = 100, these changes are 
within our errors. Similar results were obtained if each group 
was assigned the redshift of its largest component, leading 
us to conclude that the ^(u) features observed in both the 
high-redshift and low-redshift species are not related to di¬ 
vision into components, but rather reflect the underlying 
distribution of intergalactic metals. 


5 NUMERICAL SIMULATION 
5.1 Overall Properties 

To better interpret the features seen in metal absorption¬ 
line systems, we conducted a detailed smoothed particle hy¬ 
drodynamic (SPH) simulation of structure formation. Our 
goal here is to apply the same automated procedure used 
to identify metal absorbers in the LP data set to a detailed 
simulation, drawing conclusions as to what constraints our 
measurements place on the underlying distribution of IGM 
metals. For this purpose we focus our attention on a cold 
dark matter cosmological model with the same general cos¬ 


mological parameters as above, and the additional param¬ 
eters as — 0.87, and n = 1, where cg is the variance of 
linear fluctuations on the 8/i“^Mpc scale, and n is the “tilt” 
of the primordial power spectrum. The Bardeen et al. (1986) 
transfer function was used with an effective shape parame¬ 
ter of P = 0.2, and the ionizing background flux was taken 
to be a Haardt & Madau (1996) spectrum with J{v,z) = 
2.2 X 10-22 ^ ^)0.73 _ 2.3)2/i.9] ergs 

s”^ Hz”^ cm-2 sr”^, corresponding to a photoionization 
rate of 6.8 x 10"^® (1 -b 2 )° '^® exp [-(2 - 2.3)2/1.9] s’^ 

We simulated a box of size 40//i comoving Mpc on a 
side, using 320® dark matter particles and an equal num¬ 
ber of gas particles. The mass of each dark matter particle 
was 2.0 X 10® Mq, and the mass of each gas particle was 
3.4 X 10^ Mq. This yields a nominal minimum mass resolu¬ 
tion for our (dark matter) group finding of 1.0 X 10®° Mq. 
The run was started at an initial redshift of 2 = 99, and 
a fixed physical S2 softening length of 6.7 kpc was chosen, 
which is equivalent to a Plummer softening length of 2.8 
kpc. The simulation was conducted using a parallel MPI2- 
based version of the HYDRA code (Couchman et al. 1995) 
developed by the Virgo Consortium (Thacker et al. 2003). 

We used the SPH algorithm described in Thacker et al. 
(2000), although in an improvement upon earlier work, the 
maximum SPH search radius now allows us to accurately 
resolve the mean density of the box. Photoionization was 
implemented using the publicly available routines from our 
serial HYDRA code. Radiative cooling was calculated using 
the Sutherland & Dopita (1993) collisional ionization equi¬ 
librium tables, and a uniform 2% solar metallicity was as¬ 
sumed for all gas particles for cooling purposes. Integration 
to 2 = 2.0 required 9635 (unequal) steps, and four weeks of 
wall clock time on 64 processors. Outputs for post-processing 
were saved at redshifts of 2 = 8.0, 5.0, 4.0, 3.0, 2.5, and 2.0. 

From each of the final three outputs, we interpolated to 
extract two-dimensional slices of overdensity, temperature, 
and line-of-sight peculiar velocities on 24 equally spaced 
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planes. By extracting random sight-lines from each of these 
three fields, we were then able to generate simulated metal¬ 
line spectra, which could be processed in an identical fashion 
as the observed data. Before turning our attention to this 
issue, however, we first address the more basic concern of 
the overall hydrogen distribution, which serves as both a 
check of our simulation methods, and a way of fine-tuning 
the assumed ionizing background to reproduce the observed 
properties of the IGM. 


5.2 Calculation of Neutral Hydrogen Enaction 


Once the baryon density, temperature, and line-of sight ve¬ 
locity are extracted along a line of sight, constructing a sim¬ 
ulated spectrum is relatively straightforward. One obtains 
the neutral hydrogen fraction, /hi, in the IGM by solving 
the ionization equilibrium equation (Black 1981) 

a{T)npns = rci{T)nenm + J22GinHi, (9) 


where a{T) is the radiative recombination rate, Fd (T) is 
the rate of collisional ionization, J 22 is the UV background 
intensity in units of 10“^^ ergs s“^ Hz“^ cm“^ sr“^, J 22 G 1 
is the rate of photoionization, and np,ne and whi are the 
number densities of protons, electrons, and neutral hydro¬ 
gen, respectively. For the Haardt & Madau (2001) spectrum 
assumed below Gi = 2.7 x 10~^® s“^, for the original Haardt 
& Madau (1996) spectrum Gi = 3.2 x 10“'^® s“^, and for the 
(y/vm)~^ spectrum used in our simulation Gi = 3.1 x 10“^® 
s“^. For comparison, G 1 J 22 is equal to J_i 2 as defined in 
Choudhury et al. (2001) and IOG 1 J 21 as defined in Bi and 
Davidson (1997). 

If we assume that the neutral fraction of hydrogen ^ 1 
and all the helium present is in the fully-ionized form, we 
find 


/hi( a:, a) 


_ a[T(x,z)\ _ 

a [T{x,z)] + Fci [r(x, 2 )] -I- GlJ 22 ne ^( 2 ) ’ 


( 10 ) 


where the recombination rate rci(T) = 5.85 x 

exp(-157809.1/r) cm^s with T in Kelvin, 
and Black (1981) gives an approximate form for the 
recombination rate as 

a{T) _f 4.36 X (if T > 5000K) 

cm3s-i ^ I 2.17 X (if T < 5000K). 


With these expressions we can compute the neutral hy¬ 
drogen density, nm[x, 2 ( 2 :)], along a line of sight. Here x and 
2 are related by cd 2 = dx H{z), where the Hubble constant 
as a function of redshift is H{z) = + Dm(l + 2 )®. 

We choose a coordinate system such that a: = 0 at the front 
of the box and define A 2 (a;, 2 o) as the change in redshift 
from X = 0. We then construct the Lyo? optical depth as 


Taizo + Az) = 


COg 

(1 -I- 2o)yw 


/ 


nHl(x,2o) 

"" h{x,zo) 


X 


( 12 ) 


exp 


( 


^H{zo) + v{x,zo){1 + zq)—c.Az 
{ 1 +zo)b{x ,zq) 



where b(x,zo) = y (with ks the Boltzmann con¬ 

stant), nH{z) = 1.12 X 10“®(1 — y)fl6(l -I- 2 )®/i^cm“® = 
1.83x 10~’^(1-I- 2 )®cm“® (with Y the helium mass fraction), 
and ag is the Lya cross section, which can be calculated as 
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Figure 11. Top: Measured and simulated flux PDFs of the 
Lyman-a forest at two representative redshifts. Measurements 
are taken from McDonald et al. (2000) over redshift ranges of 
2.09 < 2 < 2.67 (left panel) and 2.67 < 2 < 3.39 (right 
panel), respectively, while simulations are at fixed redshifts of 
2.5 (left panel) and 3 (right panel). Bottom: Measured and simu¬ 
lated flux two-point correlation functions of the Lyman-a forest; 
/(An) =< SF{v)SF{v + Av) > where 5F = F/F—1, at two repre¬ 
sentative redshifts. The triangles are measurements by McDonald 
et al. (2000) over redshift ranges of 2.09 < 2 < 2.67 (left panel) 
and 2.67 < 2 < 3.39 (right panel) respectively and the circles 
are measurements by Croft et al. (2002) over redshift ranges of 
2.31 < 2 < 2.62 (left panel) and 2.88 < 2 < 3.25 (right panel). 
Again the simulations, represented by the solid lines, are at fixed 
redshifts of 2.5 (left panel) and 3 (right panel). 


(Tg = (37r(TT/8)^^^/Ao, (13) 

where Ao is the restframe wavelength of the transition, / is 
the appropriate oscillator strength, and ar = 6.25 x 10“^® 
is the Thomson cross section. For Lya, Ao = 1215 A and 
/ = 0.4162, which gives ag = 4.45 x 10“^®cm. With eq. II13II 
in hand, we are able to construct simulated UVES spectra 
of the Lyman-a forest by stacking vectors of optical den¬ 
sity computed from randomly extracted sightlines. These 
are then convolved with a Gaussian smoothing kernel with 
a width of 4.4 km s“^ (corresponding to the UVES reso¬ 
lution) and rebinned onto a 205,000 array of wavelengths, 
using the UVES pixelization from 3,050 to 10,430 A. Rather 
than interpolate between simulation outputs, however, we 
first turn our attention to a careful comparison between ob¬ 
servations and limited segments of spectra at fixed redshifts, 
concentrating on two main quantities: the probability distri¬ 
bution function, a single-point quantity that is sensitive to 
the overall temperature and J 22 evolution background, and 
the two-point correlation function, a measure of the spatial 
distribution of the gas. 
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5.3 Tests of The Numerical Hydrogen 
Distribution 

The probability distribution (PDF) of the transmitted flux 
was first used to study the Lyo forest by Jenkins & Ostriker 
(1991) and since then has been a widely used tool for quan¬ 
tifying the mean properties of the IGM (e.g. Rauch et al. 
1997, McDonald et al. 2000). In the upper panels of Figure 
I1 1 1 we compare the PDF as measured by McDonald et al. 
(2000) to that generated from 20 simulated spectra at rep¬ 
resentative redshifts of 2.5 and 3.0. In order to obtain the 
good agreement seen in this figure, it was necessary to ad¬ 
just the assumed J 22 values to 4.7 at « = 2.5 and 3.7 at z = 
3 (corresponding to photoionization rates of 1.3 x 10“^^ s“^ 
and 1.0 X 10“^^ s“^), down from the values of 5.4 and 4.7 
(corresponding to photoionization rates of 1.7 x 10“^^ s“^ 
and 1.5 x 10“^^ s“^) respectively, that were assumed in the 
simulations. 

This results in a slight inconsistency between the simu¬ 
lated p-T relation and the one that would have arisen if we 
had repeated the simulation with our fit values of J 22 . In 
practice however, this difference is unimportant in relation 
to our primary goal of constructing simulated metal lines. It 
is dwarfed by effects due to the uncertain evolution of the 
UV background at higher redshifts (e.g. Hui & Gnedin 1997; 
Hui & Haiman 2003), uncertainties in the normalization of 
the quasar spectra (e.g. McDonald 2000) and the extrapola¬ 
tion of the UV background from 912A to the shorter wave¬ 
lengths relevant for C iv and Si iv (e.g. Haardt & Madau 
2001). Thus our approach is more than adequate for the pur¬ 
poses of this study. With our assumed background values, 
the mean fluxes at 2 : = 2.5 and 3.0 are 0.794 and 0.692 re¬ 
spectively, while the observed values are 0.818 ± 0.012 and 
0.684 ± 0.023. 

As a second test of our simulations, we constructed the 
Lyman-a flux correlation function =< 5F{v)SF{v + 

Aw) >, which primarily provides a validation of our assumed 
primordial power spectrum P{k) and its evolution in our 
simulation. Beginning with Croft et al. (1998), the direct 
inversion of the one dimensional power-spectrum of the Lya 
flux has been seen as one of the best constraints on the shape 
of the mass power spectrum on intermediate scales (e.g. Hui 
1999; McDonald et al. 2000; Pichon et al. 2001; Groft et al. 
2002; Viel, Haehnelt, & Springel 2004). 

Again, this quantity is straightforward to extract from 
our simulated spectra. The resulting values are shown in the 
lower panels of Figure^J in which we compare them to mea¬ 
surements by McDonald et al. (2000) as well as Croft et al. 
(2002). As in the single-point case, our simulations are gen¬ 
erally in good agreement with the observed values. In fact, 
at 2 = 2.5, our simulated values are well within the range 
of values bracketed by the weakly disagreeing observational 
results. At 2 = 3.0, our simulated values provide a slight un¬ 
derprediction at small separations, although this is only just 
outside the 1-cr error in the current measurements. In sum¬ 
mary then, the gas properties of our numerical simulation 
are more than adequate to provide a Arm basis for the con¬ 
struction of Ly-a spectra, while at the same time containing 
sufficient resolution to allow us to push toward the denser 
regions associated with metal-line absorption systems. 
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Figure 12. Abundances of various species as a function of total 
hydrogen number density for a 10~^ Zq gas exposed to a Haardt 
& Madau (2001) background at 2 = 2.5 In the upper left panel, 
the temperatures corresponding to each of the curves are, from 
top to bottom, 10®'^® K (dotted), 10^ K (solid), 10^'^® K (dot- 
dashed), K (dashed), 10'^'^® K (dotted), 10® ° K (solid), 

10®.25 (dot-dashed), 10® ® K (dashed), 10®'^® K (dotted), and 
10° °K (solid). Similar labeling is used in the other panels. The 
vertical lines give the mean hydrogen number density at 2 = 3 
and 2 while the horizontal lines give the total abundance of each 
of the elements. 


6 MODELING METAL ENRICHMENT 

6.1 Calculation of Observed Metal Lines 

Extending the methods of §5.2 to construct the spatial dis¬ 
tribution of metal lines requires us to adopt an overall spec¬ 
tral shape for the ionizing background, as well as a more de¬ 
tailed calculation of the densities of various species. Here we 
assume a UV spectrum as predicted by the updated models 
of Haardt & Madau (2001; see also Haardt & Madau 1996) 
at 2 = 2.5, but shifted such that J 22 is consistent with the 
levels found in §5. Assuming local thermal equilibrium, we 
then make use of CLOUDY94 (Ferland et al. 1998; Ferland 
2000) to construct tables of the relevant species as a function 
of temperature and density at each of these redshifts, for a 
characteristic metallicity of 10' -"^0. Self-shielding in opti¬ 
cally thick regions was not taken to account. The resulting 
species densities are shown as a function of hydrogen num¬ 
ber density and temperature in Fig. II21 which is modeled 
after Figure 2 of Rauch, Haehenelt, & Steinmetz (1997). 

In this figure, we see that, roughly speaking, C iv traces 
the widest range of environments, while Si iv. Mg ii, and 
Fe II probe progressively denser regions. Thus while an ap- 
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Figure 13. Correlation functions and column density distributions for models in which metallicity is assumed to be constant throughout 
the simulation, or a simple function of density. The filled points in the upper panels show high-metallicity models in which Z = 10~^Zq 
and Z = while the filled points in the lower panels show lower metallicity models in which Z = 10~^Zq and Z = 

Fifty-seven simulated QSO sightlines were averaged in the high-metallicity models and 114 were averaged in the lower- 
metallicity cases. The open circles give the measured C iv correlation function, and the dashed lines give the fit to the column density 
distribution as in Figure Q 


preciable level of C iv is found in only a few times overdense 
2 = 3 gas, comparable levels of Si iv are achieved only in 
denser regions with A = p{x)lp « 10; and while Mg ii is 
found at similar densities to Si iv, Fe ll is only detectable 
in A > 100 regions, orders of magnitude denser than most 
C IV regions. 

Similarly, C iv is detectable over a large temperature 
range, covering from 10^ K up to « 10® K. While Si iv is 
also relatively stable with respect to temperature changes. 
Mg II and Fe ii are much more fragile, and their abundances 
fall away quickly above ~ 10® K. From these results, we see 
that the correct modeling of Si iv requires simulations that 
probe to densities 10 times higher than those most rele¬ 
vant to C IV, although 5civ(u) and 5siiv(u) trace each other 
closely. Thus while our numerical modeling was carried out 
at the highest resolution possible, we nevertheless limit our 
comparisons to C iv to minimize any remaining numerical 
effects. 


6.2 A Nonlocal Dependence 

Having determined the number densities of each of the 
species of interest as a function of temperature and density 
in a 10-2^0 medium, we then applied these calculations to 
extract simulated metal absorption spectra from our simu¬ 


lations. As a first step, following Rauch, Haehnelt, & Stein- 
metz (1997), we assumed a constant metallicity across the 
simulation volume and extracted sightlines of rc iv using an 
appropriately modified version of eq. m-- 


TCiv,i{zo + A2) = 


CO'CiVj'i f ^ iv (x,zo) 
(l-l- 2 o)v^y 6civ(a:,2o) 


(14) 


exp 


( xH(zq)+v(x,zq)(1 + zq) — cAz 

(l+ZQ)bc iv (a:,zo) 



where now &civ(a:, 20 )= y —nciv( 2 ) is the mean 
C IV density, and (Tciv.i is the cross section corresponding 
to the ith absorption line of the C iv doublet. These we 
compute directly from eq. unii taking (Ao,civ,i, Ao,civ, 2 ) — 
(1548.2,1550.8) and (/civ.i,/civ. 2 ) = (0.1908,0.09522). For 
the low metallicities relevant for the IGM, the effects of 
changing metallicity can be modeled as a simple linear shift 
in the species under consideration. 

In contrast with the fixed-redshift comparisons de¬ 
scribed in §5, our goal was to construct simulated data sets 
that corresponded as closely as possible to the full ESO 
Large Program data set. In this case, instead of stacking 
together spectra drawn from a single output, we instead al¬ 
lowed for redshift evolution: drawing slices from the output 
that most closely corresponded to the redshift in question. 
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taking n oc (I + 2 )®, and interpolating between CLOUDY ta¬ 
bles with appropriate J 22 values. Finally, we applied Poisson 
noise corresponding to a signal-to-noise of 100 per pixel. 

Each spectrum generated by this method was subject to 
the same two-step identification procedure that was applied 
to the real data, and the resulting fits were subject to the 
same N and b cuts described in §2.2. The line lists compiled 
in this way were then used to generate correlation functions 
and column density distributions that directly parallel those 
calculated from the Large Program data set. 

These are shown in Fig. 13, in which we explore a low- 
metallicity model (10’ -^Zq) roughly consistent with pre¬ 
vious estimates (e.g. Rauch, Haehnelt, & Steinmetz 1997; 
Schaye et al. 2003) as well a higher-metallicity model 
(10“^Z©). Note that at these metallicities, changes in Z 
can be modeled simply by boosting the C iv density derived 
from the CLOUDY tables by a linear factor. Increasing the 
metallicity by a factor of 10 in this way has very little ef¬ 
fect on the correlation function: decreasing ^civ(w) in the 20 
and 35 km s“^ bins by roughly a factor of 1.5, while leaving 
the rest of the correlation function largely unchanged. In all 
cases these values fall far short of the clustering levels seen in 
our observational data, and they lack the conspicuous bend 
observed at « 150 km s“^. 

However, changing metallicity has a large effect on 
the column density distribution. The low-metallicity model 
is consistent with observations over the range of 12.5 < 
log Y[cm“^] < 13.5, and slightly overpredicts the number of 
large components (which have a negligible impact on the cor¬ 
relation function). The high-metallicity model, on the other 
hand, overpredicts the number of components for all column 
densities log Y[cm“^] > 13.0. 

The poor fit to the correlation function is perhaps not 
surprising given the known inhomogeneity of the IGM metal 
distribution (e.g. Rauch, Haehnelt, & Steinmetz 1997). Most 
recently, this has been quantified by Schaye et al. (2003), 
who applied a pixel optical depth method to derive a non¬ 
linear relation between the local overdensity A of hydrogen 
and the local carbon abundance. Over a large range of envi¬ 
ronments, they found [C/H] oc with a large variance. 
Is it possible that accounting for this relationship would be 
able to resolve the discrepancy in iv(«)? In order to ad¬ 
dress this question, we repeated our experiment, assuming 
that the local density followed the best-fit relationship de¬ 
rived by Schaye et al. (2003). Again we considered both high 
and low metallicity models, resulting in correlation functions 
and column density distributions that are shown in Fig. 13. 

As our results, which depend on a line-fitting procedure, 
are biased to the densest regions, the “zero point” metallic¬ 
ity of these models are naturally shifted to lower values. 
Thus, the A^^^IO'^^© and the Zq models shown 

in these figures yield similar numbers of components as the 
single metallicity 10“^^© and 10~^Zq models, respectively. 
In particular, the lower metallicity model allows us to obtain 
good agreement with the observed column density distribu¬ 
tion, while assuming mean metallicity values more in line 
with previous estimates (e.g. Hellsten et al. 1997; Rauch, 
Haehenelt, & Steinmetz 1997; Schaye et al. 2003). 

Introducing a A dependence has almost no effect, how¬ 
ever, on the correlation function, neither boosting it at low 
separations nor introducing a feature at « 150 km s“^. Thus 
is appears that this nonlinear relationship is not the source of 


the clustering properties of the metal-line components, and 
rather that the large variance seen in the pixel-by-pixel re¬ 
sults hides a third parameter that determines these features. 
In fact, in Paper I, we described just such a key parameter, 
the separation from a large dark-matter halo. 

7 SOURCES OF INTERGALACTIC CARBON 

IV 

7.1 Distribution of Metals and Identification of 
Sources 

While the observed features in the C iv correlation function 
can not be understood in terms of a local nonlinear rela¬ 
tionship between the metal and density distributions, we 
saw in Paper I that these features could be easily explained 
in terms of the distribution of metal sources. In that work 
we used a pure dark-matter model to describe C iv compo¬ 
nents as contained within bubbles centered around sources, 
and we interpreted the amplitude and the knee in the C iv 
correlation function in terms of the source mass and bubble 
size, respectively. In this investigation we develop a similar 
model, but make use of the full gas and CLOUDY modeling 
described in §5 and §6. 

Following Paper I, we adopt a parameterization in 
which all metals are found within a comoving radius Ra of 
a dark matter halo whose mass is above a fixed value, Ms. 
To facilitate comparison with our previous modeling, as well 
as to allow for future comparisons with analytic approaches, 
we identified all sources at a fixed redshift of « = 3. In par¬ 
ticular, halo detection was performed using the HOP algo¬ 
rithm (Eisenstein & Hut 1998) with parameters dpeak = 160, 
(^saddle = 140, and Jouter = 80. The centres of these groups 
were then traced forward in time to the z = 2.5 and 2.0 
slices such that exactly the same groups could be selected 
from all the simulation slices, accounting for appropriate pe¬ 
culiar motions. 

As in Paper I, our choice of 2 : = 3 is not meant to imply 
that enrichment occurred at this redshift, but rather that it 
occurred at an unknown redshift higher than the observed 
range, centered on groups whose large-scale clustering was 
equivalent to 2 = 3 objects of mass Ms. For each choice 
of Rs and Ms, we then painted bubbles of metals on our 
simulations, as illustrated in Fig. 14. While increasing Ra 
has the obvious effect of increasing the volume filling factor 
of metals, increasing Ms not only lowers this filling factor, 
but also clusters the bubbles more strongly. This can be 
most easily seen by comparing models with similar filling 
factors. For example, comparing the Ms = 1 x 10^^ M©, 
Ra = 1.6 comoving Mpc slice to the Ms = 5x 10^^ M©, Ra = 
2.4 comoving Mpc slice indicates that a similar fraction is 
enriched with metals in both cases, but these regions are 
spread over a considerable area in the lower-mass case and 
concentrated into dense knots in the higher-mass case. 

7.2 Properties of Carbon-IV 

From slices such as those shown in Fig. 14 we were able 
to generate simulated QSO absorption spectra, in a manner 
exactly parallel to that described in §6.2: drawing lines-of- 
sight for the various time outputs, piecing them together by 
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Figure 14. The C iv distribution of a ^ = 3 slice in the simulation, dark regions are those contained within a distance Rs from a dark 
matter halo of mass Ms, with parameters as labeled in the panels. These regions are used in constructing simulated spectra, while all 
gas outside them is considered to be metal-free. From left to right, and then top to bottom, the overall volume filling factors of these 
models are: 3.3%, 8.6%, 16.1%; 3.4%, 8.6%, 16.7%; and 1.7%, 5.5%, 11.6% 


evolving the mean density, interpolating between CLOUDY 
tables, and applying realistic levels of Poisson noise. In this 
case the metallicity was assumed to be at a fixed value Zf, 
within each bubble, and zero everywhere else. Note, however, 
that our measurements are insensitive to C iv components 
with columns below 10^^ cm“^, a thus a more widely dis¬ 
persed, lower-level contribution to IGM metals (e.g. Schaye 
et al. 2003; Bergeron & Herbert-Fort 2005) can not be ex¬ 
cluded. 

In Paper I, our modeling made use of a parameter b 
that controlled the impact parameter associated with each 
sub-clump within a bubble. In our more physical modeling, 
this role is played by Zb, which we fixed at an initial value 
of 1/2O.^0. These spectra were again analyzed by our au¬ 
tomated procedure, and in Figures 15 and 16 we compare 
the resulting correlation functions and column density distri¬ 
butions with those measured in the Large-Programme data 
set. 

These plots reflect the trends seen in the slices. In¬ 
creasing the mass concentrates the metal into fewer regions, 
boosting the correlation function, particularly at large sep¬ 
arations. Increasing Rs, on the other hand, impacts the cor¬ 
relation function primarily at smaller separations, and has a 


strong impact on the total number of C iv components de¬ 
tected per spectrum. Prom Fig. 15, the best-fit models are 
the Ms = 5x 10^^ Mq k Rs = 2.4 Mpc, Ms = 10^^ Mq & 
Rs = 2.4 Mpc, and Ms = 10^^ M© & = 3.2 Mpc cases, 

with filling factors of 5.5%, 8.6%, and 11.6% respectively 
although several of the lower filling-factor cases produce so 
few lines as to be difficult to evaluate in detail. Similar fill¬ 
ing factors have been advocated by Fieri & Haehnelt (2005) 
on the basis of O vi measurements. The large Ms values 
we derive are also suggestive of the regions around Lyman- 
Break galaxies (LBGs), which are observed to be clustered 
like « 10^^ ® Mq halos at « = 3 (Porciani & Giavalisco 2002), 
and for which a strong cross-correlation with C iv absorbers 
has been measured (Adelberger et al. 2003). It is also rem¬ 
iniscent of the association between galaxies and C iv ab¬ 
sorbers put forward in BSR03. 

In Fig. 16 we better quantify the number of components 
in each model, by constructing simulated column density 
distributions as discussed in §3.1. Here we see that regardless 
of our choice of source mass and bubble radius, all these 
models fall short of reproducing the observations. Due to 
the relatively small filling factors of such bubble models, our 
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Figure 15. Correlation function of simulated C iv components, with an assumed bubble metallicity of 1/20Zq. Panels are labeled by 
their assumed Ms in units of Mq and Rs in units of comoving Mpc. In each panel the open circles summarize the observational results, 
while the filled circles represent the experimental results, as averaged over 114 spectra. Each panel is labeled by the average number of 
C IV components detected per simulated spectrum. 


choice of Zi, = 1/20Zq is not able to generate the relatively 
large number of C iv absorbers seen in the data. 

In order to improve this agreement then, we considered 
a model in which we assume a higher bubble metallicity of 
Z), ~ 0.2Zq, generating the 5civ('u) and log/(Ai") values seen 
in Figures 17 and 18. As was the case for our b parameter 
in Paper I, varying Zi, has relatively little impact on the 
C IV correlation function, although the increased number of 
components does result in less noise. 

Thus the high-metallicity simulations display the same 
trends and best-fit models as were seen in the lower Zb case, 
however the improved signal allows us to distinguish the 
Ms = 10^^ Mq, R = 2.4 Mpc model as a somewhat better 
match to the data than the Ms = 5 x 10^^ Mq, R — 1.6 Mpc 
and Ms = 10^^ Mq, R = 3.2 Mpc models. The improved sig¬ 
nal in Fig. 17 also enables us to reject cases with very low 
filling factors. In particular, we see that the models with the 
smallest bubble sizes do not reproduce the observed 150 
km s“^ elbow, exceeding the measured ^cw{v) at small sep¬ 
arations. Furthermore, models with Ms = 10*^^ M© are now 
seen to fall far short of the observed correlation function at 
large separations, particularly if we consider the models with 
Ra > 0.8, which are not overly peaked at small distances. 

Finally, assuming a mean bubble metallicity of 1/5^0 
has a large impact on the column density distribution, ap¬ 
proximately doubling the number of detected components 
and bringing our best-fit model in rough agreement with 
the data, although perhaps even this value is slightly low 


in our best-fit cases. It is clear that we are forced toward 
these values because much of the gas around « IO'^^Mq is 
heated by infall above « 10®'® K, and thus it is largely in the 
outskirts of our bubbles in which C iv absorbers are found. 


While, at face value, this metallicity is widely discrepant 
with other estimates, there are nevertheless two reasons to 
take it seriously. Firstly, the dense metal-rich regions in 
our model are observed to be enriched to similar levels at 
z; = 1.2. At this point, the LBG-scale halos around which 
we have placed our metals are expected to have fallen into 
clusters, and thus the IGM gas is detectable through X-ray 
emission in the intracluster medium (ICM). In fact, detailed 
Chandra and XMM-Newton observations indicate ICM iron 
levels oi Z — 0.20to.o5 at 2 = 1.2 (Tozzi et al. 2003), im¬ 
plying that at even higher redshifts these metals have been 
efficiently mixed into the diffuse gas that forms into clus¬ 
ters. Secondly, we note that more than enough star forma¬ 
tion has occurred by 2 = 2.3 to enrich these regions to our 
assumed values. Indeed, comparisons between the integrated 
star formation history and more standard estimates of IGM 
metallicity have shown that a large fraction of 2 « 2 metals 
have so-far escaped detection (Pettini 1999). Thus, we find 
no compelling reason to dismiss this high metallicity value 
as spurious, although we emphasize that it has no direct 
impact on our derived clustering masses and bubble sizes. 
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Figure 16. Column density distribution of simulated C iv absorption components, with an assumed bubble metallicity of 1/20^0. 
Models are as in Fig. 16 and in each panel the filled points represent the simulation results, while the dashed line is the fit given in Figure 
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8 AN ANALYTIC MODEL 

While our simulated bubble model provides a compelling 
picture of the C iv and Si iv distribution at z « 2 — 3, 
it leaves open the question as to properties of Mg ii and 
Fe II. Yet the detailed modeling of these species is beyond 
the capabilities of our simulation. As we saw in §6.1, the en¬ 
vironments of Mg II and Fe ii are denser than C iv and Si iv, 
particularly in the case of Fe ll. Even more constraining is 
the fact that almost all our detections of these systems fall 
well below our final simulated redshift of 2, with the major¬ 
ity lying in the 0.5 ^ z ^ 1.5 redshift range. 

In Fig. 1^ we saw that while the overall shape and cor¬ 
relation function of these species is comparable to that ob¬ 
served in C IV and Si iv, the magnitude and long-separation 
tail of 5(u) are shifted upwards by a factor associated with 
the cosmological growth of structure. While reproducing 
these trends is beyond the capabilities of our simulation, 
they can nevertheless be studied from an approximate ana¬ 
lytic perspective. 


8.1 Derivation 

In §7 we found good agreement between our observations 
and a model in which we painted metals around biased re¬ 
gions in our simulation. Analytically this corresponds to 
a picture in which the metal lines observed at Zobs come 
from clumps that are within a fixed radius of the sources 
of IGM metals. These pollution sources are associated with 


relatively rare objects of mass Mp that ejected metals into 
surroundings at a high redshift Zp > Zobs. After enrichment 
these components continue to cluster gravitationally to Zobs. 

In the numerical simulations, these pollution centres are 
identified with the galaxies of mass Ms at a redshift of z = 3. 
However, this mass and redshift were intended only to quan¬ 
tify the bias of sources, and it is more probable that they are 
really related to less massive, higher redshift objects, which 
exhibit similar clustering properties (Scannapieco 2005). 

Let us consider then four points: the centres of two 
clumps (1 and 2), which we observe as metal-line compo¬ 
nents, and the centres of two bubbles (3 and 4), which cor¬ 
respond to the sources of pollution. We require that the pol¬ 
lution sources correspond to peaks \i.e. linear overdensities 
with a contrast larger than 5cv = 1 . 6811 ( 2 ;^)“^] at a redshift 
Zp > Zobs and at a mass scale Mp. The clumps, on the other 
hand, are associated with the C iv absorbers themselves, 
and correspond to peaks at a mass scale Me. In the linear 
approximation, these fields satisfy a joint Gaussian proba¬ 
bility distribution which is specified by the block correlation 
matrix: 


M 


Ccc(o) ^cc(fl 2 ) 

Ccc(l’l2) Ccc(O) 
Cep (»'13) Cep (^■23) 

_ Cep (ri4) Cep (r24) 

]^cc Cep 

Cpc Mpp ’ 


Cep (ris) Cep (ri4) 1 

Cep (r23) Cep (r24) 

Cpp(O) Cpp(^34) 
Cpp(^34) Cpp(O) . 


(15) 


where Vij = ||ri — rjH, and ^pp, ^cc and ^cp refer respec- 
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Figure 17. Correlation function of simulated C iv components, with an assumed bubble metallicity of 1/5Zq. Panels and symbols are 
as in Fig. 16. 


tively to the correlation between pollution centres, satellite 
clumps, and the cross correlation between clumps and the 
centres. The joint probability of having a peak of an ampli¬ 
tude larger than Sct at the four points is given by 


p(l,2,3,4) = 


X 


X 


1 

dyr^y^det |M| 




d54 


_ (<51,^2'^3’*^ 4 )^ 

e 2 


(16) 


We will evaluate this expression, assuming that the thresh¬ 
old that defines the object is high relative to the correspond¬ 
ing rms densities and taking the correlation between the 
metal line clumps and the centres of pollution to be small. 
We shall not assume the smallness of the centre-centre nor 
clump-clump correlation, of which the first is the most im¬ 
portant. In this limit 


det |M| « 

and 

p(l,2,3,4) 


Mcc Mcc ^cpMpp 

-MppCpcM”! Mpp 

det |Mcc| • det |Mpp| , 


(17) 

(18) 


47r2 det | Mcc I \/det | Mpp | 


d^l / ( 1^26 

Js^r 


X 


X 


d53 / d54e-5(^3,^4)- M-b(^3,^,) 


0(^1 ,^ 2 CepMpp^ .(153.154) 


(19) 


In the high peak limit, the last cross-correlation term can 
be factored out from the integrals (see Appendix) yielding 

p(l, 2, 3, 4) Ri p(l, 2)p(3, ( 20 ) 


where p(l, 2) and p(3, 4) are computed from eq. ifm . Or, 
explicitly. 


p(l,2,3,4) « 7^i^cc^t'pp^C[c(ri2),J^cc]C[c(r34),^'pp] (21) 


47J-2 

X exp 


+ z^cc^p 


X 


1 + Ccc(ri 2 ) 1 + Cpp(r34) 

Ccp(y'l3) + Ccp(y'24) + Ccp(?"l4) + Cep S 23 ) 
[1 + ‘^PP ( 7 - 34 )] [1 -b Ccc(ri2)] 


where the function C{x, v) is defined in the Ap¬ 
pendix, and we define the cross-correlation coefficients 

as Ccc = Ccc(r)/^cc( 0 ), Cpp = Cpp(^)/Cpp( 0 )i Cep = 

^cp(c)/-\/Ccc(0)Cpp(0), and the normalized density thresh¬ 
olds as P'ce = ^cr j \l^cc (0) and ^Zpp = dcr/ y/^pp (0) * 


^ The cross-correlation coefficients ^cc (r)/?cc (0) and 
^pp(x)/^pp(0) reach unity at r = 0 and thus cannot be assumed 
small everywhere. At the same time icp{r)/\J ?cc(0)^pp(0) is 
always less than unity if Me and Mp do not coincide (Schwartz 
Inequality). In particular the smaller its maximum value. 
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Figure 18. Column density distribution of simulated C iv absorption components, with an assumed bubble metallicity of l/SZ©. Panels 
and symbols are as in Fig. 17. 


Note that in eq. 121 ll the correlation functions in the 
denominator are not assumed to be small, which allows for 
proper accounting of the case when two clumps or two pollu¬ 
tion centres are the same. For example, setting ri 2 = r 34 = 0 
properly recovers the bivariant joint probability p(l,3) to 
find a clump at a separation ri 3 from the centre of pollution 
(equal in this case to ri 4 = r 23 = r24). 

In our model only those clumps that lie within the 
spherical bubble around some pollution centre are observed 
to have metals. The correlation function of clumps of mass 
Me that are within spherical bubbles around peaks corre¬ 
sponding to the mass Mp, is defined as 

t(<5i, (53)j3((52, < 54 ) [l-b f(ri2)] = p((5i, (52, ^3, < 54 ), (22) 

where the bar denotes averaging over position of pollution 
centres within a distance Rs around two metal-rich clumps 
at a fixed separation ri 2 . Note that our definition of the cor¬ 
relation function, C(ri 2 ), is not equivalent to the estimator of 
the underlying correlation function of all the clumps of mass 
Me, ^(ri 2 ), nor is it equivalent to the conditional correlation 
function if there were a source of metals (a high peak of the 
scale Mp) in the vicinity of every small halo, ^e(ri 2 ), which 
would be given by 

p(5i)p((52)p((53,54) [ 1 -b ^c(ri2)] =p((5i, 52,^3, 54 ). (23) 

Furthermore, ^(’’ 12 ) depends on the underlying two-point 

achieved at r = 0, the larger the difference between the scales 
describing the clumps and the pollution centres. 


correlation of small clumps, the correlation of the sources, 
and the cross-correlation between clumps and sources. This 
last term is subject to the most modification should the 
physics of metal dispersal change. However it mostly affects 
the biased density of small clumps in the vicinity of the 
sources relative to the cosmological mean, which is precisely 
the excess factored out in eq. 12211. 

Thus eq. 12211 describes the correlation of metal compo¬ 
nents at the redshift of pollution, which is dominated by the 
clustering of the pollution sources. Subsequent gravitational 
clustering of enriched metals then leads to further amplifi¬ 
cation of the correlation in the linear approximation as 

^(r,Zobs) = l-D(zobs)/-D(z:p)f ^(r,Zp), (24) 

where D(z) is the linear density growth factor. This growth 
is suggestive of the difference between the C iv and Mg ii 
correlation functions, as we saw in FigEl as well as the hints 
of evolution seen in 5civ(w) and 5Mgii(u) in Figs. |21and|3 


8.2 Application to observed metal absorbers 

In Figure Eni we fit our analytic model to the data. In the 
left panel we adopt the parameters used in our numer¬ 
ical simulations, identifying metal pollution centres with 
Mp — 10^^ Mq objects at a redshift Zp = 3 and metal 
rich clumps with collapsed objects of Me = 10®M©, with 
Rs — 2.4 comoving Mpc. At Zobs = 2.3, the analytic fit re¬ 
produces the measured Cciv(u) at large velocity separations, 
where it is dominated by the correlation between pollution 
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Figure 19. Comparison of our analytic model with the data. Left 
panel: Low-redshift model for metal sources. Dotted lines repre¬ 
sent the linear clustering at Zobs = 2.3 (lower) and Zobs = 1-15 
(upper) of clumps observed in the vicinity of the pollution cen¬ 
tres, with Mp = Me = 10®Mq and Zp = 3. Solid lines 

show the effect of applying a non-linear correction to these mod¬ 
els. Finally, the dashed curve shows a nonlinear Zobs = 1-15 model 
in which Mp = IO^^Mq, and Zp = 3, but now Me = 

Right panel: Lower solid curve corresponds to Mp = 3 x 10® M©, 
Me = 1O®M0, Zpoi = 7.5, Zobs = 2.3, with no nonlinear correc¬ 
tion applied. The upper solid line is a further linear extrapolation 
of this model to Zobs = 1.15. The upper dashed line is a linear 
model again with Zobs = 1-15 Mp = 3 x 10®Mq, Zpoi = 7.5, but 
with Me raised to 10® Mq. For all curves the comoving size of the 
bubble is 2.4 Mpc. 


centres, but it falls short at small velocities, where ^civ(u) 
is dominated by the clump distribution within each bubble. 

This is because the smoothing imposed by choosing 
Mp « 10^^Mq is similar to the 2.4 Mpc bubble size, and 
thus our linear formalism is insufficient to describe distances 
less than Rs- In reality, the nonlinear collapse of Mp would 
have moved in new material to fill in this region. To mimic 
such nonlinear effects at small radii, we apply the prescrip¬ 
tion 5cr —> (5cr + (1 — l/5cr)^ (Mo & White 1996), resulting 
in the dashed curve. This correction, while crude, is seen to 
recover a ^civ(u) that is similar to our simulated 10^®Mq, 
Rs = 2.4 Mpc case (and thus the observed correlation func¬ 
tion), confirming that the discrepancy at small distances is 
caused by our neglect of nonlinear motions. 

Next we turn our attention to Mg ll and Fe ll, which are 
observed at lower redshifts z « 1.2. As we saw in Figure El 
the rise of the correlation amplitude of these species relative 
to C IV and Si iv is generally in agreement with the hypoth¬ 
esis of linear growth of gravitational clustering of a fixed 
population of objects from z = 2.3 to z = 1.15, although 
there are significant discrepancies at small radii. Again we 


plot both a linear Zpoi = 3, Mp = 10^®M©, Me = 10®M© 
model observed at Zoba = 1.15 and a similar model in which 
a nonlinear correction has been applied. While the nonlinear 
curve does well at most radii, a shortfall is seen at z 100 
km s“^, similar to the discrepancy between the “shifted” 
Cciv(w) curve and the 5Mgii(r’) curves in Figure El Based 
on our plots of the species fraction as a function of envi¬ 
ronment, an important difference between these species is 
clear. As Mg ii can only survive in regions with a low ion¬ 
ization parameter, it is biased toward much denser regions 
than C IV, which correspond in our analytic models to higher 
clump masses. Raising Me to 10^®M© to account for this ef¬ 
fect leads to the dashed curve in the left panel, which again 
agrees well with the data. 

As discussed above, however, it is likely that the origin 
of metal pollution lies at higher redshift from sources of a 
lower mass, whose comoving clustering properties are iden¬ 
tical to M « 10^®M© galaxies identified at z = 3. Indeed, 
these biased high-redshift sources may be the progenitors 
that later grew into large z = 3 galaxies. In the right panel 
of Figure^lwe explore such a high redshift model, in which 
we take Zpoi = 7.5 and Mp = 3 x 10®M©, so that the bias 
of our sources is the same as for Mp — 10^® M©, objects 
forming at Zp = 3. Adopting a similar Ve as in the Zp = 3 
case results in the solid curves. As the smoothing due to the 
Lagrangian radius associated with Mp = 3 x 10® is minimal, 
no nonlinear correction is necessary and our simple model 
provides a reasonable fit to the C iv and Si iv components 
observed at z = 2.3. 

Similarly, extrapolation of the same objects to z = 1.15, 
the mean redshift for Mg ii and Fe ii systems, matches their 
large scale (v > 300 km s“^) correlations quite well, al¬ 
though the data at these separations is sparse. At small ve¬ 
locities the difference between the environments of the two 
species becomes important, and linear scaling does not com¬ 
pletely explain the enhancement of correlation amplitude 
in Fe ii and Mg ii relative to C iv and Si iv. As in the 
low redshift case, if we associate these species with larger 
clumps, the fit is much improved at small radii, resulting in 
the dashed curve. 

In summary, our simple analytic model generally repro¬ 
duces the features seen in our simulations of the C iv and 
Si IV correlation functions, although a nonlinear correction 
is necessary in the Zp = 3 model. Linearly extrapolating 
these models to lower redshift results in a good fit to Mg ll 
and Fe ii at large distances, although a fit at smaller dis¬ 
tances requires us to use larger clump masses, associated 
with denser environments. Finally, we find that there is a 
strong degeneracy between Mp and Zpoi, with a family of 
sources with similar biases producing acceptable fits. 


9 CONCLUSIONS 

While intergalactic metals are ubiquitous, the details of how 
these elements made their way into the most tenuous regions 
of space remains unknown. In this study we have used a 
uniquely large, homogeneous, and high signal-to noise sam¬ 
ple of QSO sightlines to pin down the spatial distribution of 
these metals and combined this with advanced automated 
detection techniques and a high-resolution SPH simulation 
to pin down just what we can learn from this distribution. 
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Our study has been focused on four key species: C iv and 
Si IV, which serve as tracers of somewhat-overdense regions 
from redshifts 1.5 to 3.0, and Mg ii and Fe ii, which trace 
dense, lower-redshift (z = 0.5 — 2.0) environments. No evo¬ 
lution in the column density distributions of any of these 
species is detected. 

In the high-redshift case, C iv and Si iv trace each other 
closely. For both species, ^(v), exhibits a steep decline at 
large separations, which is roughly consistent with the slope 
of the ACDM matter correlation function and the spatial 
clustering of a « 3 Lyman break galaxies. At separations 
below « 150 km s“^, this function flattens out considerably, 
reaching a value of ^(v) « 50, if u 50 km s“^. Our data 
also suggests that 5civ('y) evolves weakly with redshift, at a 
level consistent with the linear growth of structure. 

The distribution of metals as traced by ^civ(u) is ex¬ 
tremely robust. We And that it remains almost completely 
unchanged when minimum or maximum column density cuts 
are applied to our sample, even if they are so extreme as 
to eliminate over two-thirds of the components. We have 
also linked together C iv components into systems, using 
a one-dimensional friends-of-friends algorithm, with linking 
lengths of uimk = 25, 50, and 100 km s“^. In all cases, 
the line-of-sight correlation function of the resulting systems 
matches the original component correlation function (within 
measurements errors) at separations above uunk. Finally, the 
Si IV/C IV ratio shows no clear dependence when binned as 
a function of separation, suggesting that the features seen 
in ^civ(v) and Csiiv(w) do not result from fluctuations in the 
ionizing background. 

Thus none of our tests indicate that the observed dis¬ 
tributions of C IV and Si iv represent anything but the dis¬ 
tribution of intergalactic metals at z = 1.5 — 3.0. This moti¬ 
vated us to carry out a confrontation between our C iv obser¬ 
vations and detailed simulations of IGM metal enrichment, 
which paralleled previous comparisons for the Lyman-alpha 
forest. Furthermore, the advanced automatic-detection pro¬ 
cedures described in §2.2 (see also Aracil et al. 2005) allowed 
us not only to compare simulated and observed spectra, but 
generate simulated line lists in a manner that exactly paral¬ 
leled the observations. 

Using these tools, we found that the observed features 
of the C IV line-of-sight correlation function can not be re¬ 
produced if the IGM metallicity is constant. Rather any such 
model falls far short of the observed iv (u) amplitude and 
fails to reproduce flattening seen below « 150 km s“^. Fur¬ 
thermore, adopting a local relation between overdensity and 
metallicity, as observed by Schaye et al. (2003), has little or 
no effect on these results. 

On the other hand, rough agreement between simulated 
and observed C iv correlations is obtained in a model in 
which only a fraction of the IGM is enriched. Emulating 
the simple model in Paper I, we explored a range of models 
in which metals were confined within bubbles of radius Ra 
about 2 = 3 sources of mass Ms, where these quantities are 
not meant literally as source redshifts and masses, but rather 
as tracers of the bias of the 2 poi > 3 source population. 
Varying these quantities, we derived parameters that suggest 
large metal bubbles, Rs ~ 2 comoving Mpc, around highly- 
biased sources, with Ma ~ 10^^ M©. 

These results are suggestive of the association between 
galaxies and G iv absorbers put forward in BSR03, and the 


high cross-correlation between LBGs and G iv absorbers 
measured by Adelberger et al. (2003). Yet this does not mean 
that LBGs are the sources of intergalactic metals, only that 
the 2poi > 3 sources were biased like LBGs. In fact, the 
case for outflows escaping lower-redshift starbursts is much 
more convincing than for dwarves (Martin 2005). Similarly, 
Ra need not be interpreted as the ejection radius of each 
source, but instead as the distance at which bubbles from 
multiple sources overlap. Our best fit Ma and Ra values are 
independent of the assumed bubble metallicity, although the 
low ^ 10% volume Ailing factors of these models forces us 
to use large « ^/^Zq values to reproduce the observed C iv 
column density distribution. Note however that given the 
high bias of our enriched regions, such metallicities may be 
necessary to reconcile 2 « 2.3 measurements with 2 « 1.2 
observations of the iron content of the IGM in high-redshift 
galaxy clusters (Tozzi et al. 2003). 

At lower redshifts, the line-of-sight correlation functions 
of Mg II and Fe ii are consistent with the same enriched 
regions seen in G iv and Si iv, but now “passively” evolved 
down to 2 « 1.2. Again both ^Mgii(u) and ^Feii(v) trace 
each other closely, and exhibit the same steep decline at 
large separations and flattening at small separations as were 
seen in ^civ(u) and ^siiv(v). Also as in the higher redshift 
case, the Mg ii correlation function remains unchanged when 
minimum and maximum column density cuts are applied, 
and linking together Mg ii components into systems has no 
strong impact on ^Mgii(u) outside separations corresponding 
to the linking length. 

Although Mg II and Fe ll are detected in regions that 
can not be simulated numerically, we are nevertheless able 
to develop an analytic model that allows for a simple anal¬ 
ysis of these species. Testing our model against ^civ{v) and 
Csiiv(f), we And generally good agreement with the data 
for similar values of mass and Ra as in the numerical case. 
Pushing the model to lower redshift, we And that the same 
parameters do well at reproducing the clustering proper¬ 
ties of Mg II and Fe ll, especially when we account for the 
fact that the species are found in denser environments. Fi¬ 
nally, we also And agreement with the observed ^civ(u) and 
Csiiv('y) at 2 = 2.3 and ^Mgii('y) and ^Feii(i’) at 2 = 1.15 
and a high-redshift analytic model in which Zp = 7.5 and 
Mp = 3 X 10® M©, illustrating the strong degeneracy be¬ 
tween Mp and Zp for similarly biased sources. 

Taken together, our 2 « 2.3 and 2 « 1.2 measurements, 
numerical simulations, and analytic modeling paint a con¬ 
sistent picture of IGM enrichment. The distribution of inter¬ 
galactic metals does not appear uniform, nor simply depen¬ 
dent on the local density, but rather it bears the signature 
of the population from which it came. While the 2 > 3 red¬ 
shift of metal ejection is unknown, a joint constraint on the 
masses and redshifts of the objects responsible for IGM pol¬ 
lution remains compelling. Models of IGM enrichment must 
come to terms with the observed biased sources of inter¬ 
galactic metals. 
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APPENDIX A: TWO POINT JOINT PROBABILITY DISTRIBUTION OF HIGH PEAKS IN 
GAUSSIAN FIELDS 

In this section we present formulas for the joint two point probability for the peaks in the Gaussian field that are used in our 
analytic model. These allow for the pair of peaks to have different scales and improve on the asymptotic results for the high 
peaks. 

In the standard cosmological picture one identifies a collapsed object of mass M with a peak of height 5 > Scr in the 
density field 5{x), smoothed with a top-hat window function W{R) with the scale R = In the limit of large height 

the geometrical peaks of the Gaussian field can be approximately described as just the regions of high field values. This is the 
approximation that we adopt. 

We shall need, first, the variance of the smoothed density field 


a^= P{k)W^{kR)k‘^dk, 


(Al) 


where P{k) is the power spectrum of the density field and the Fourier image of the top-hat window is 
sin kR cos kR 


W{kR) = 47ri^'' 


(A2) 


(Ai?)3 (KR)^ 

and, second, the correlation function between the values of the field at two positions, separated by the distance ri 2 = xi — X 2 


^(ri2) = / P{k)^lTj^^W{kRi)W{kR2)k^dk, 
' kri2 


(A3) 


where the value at point 1 is taken after the field is smoothed on a scale Ri while at the point 2 the field is evaluated after 
smoothing on a scale i? 2 . If Ri = R 2 , then ^( 0 ) = while in general ^( 0 ) < cricr 2 . 

To evaluate the probability distribution functions used in §8 we begin with the well-known result for the one point 
probability of the field height to exceed &cr ■ 


p(l) = 


y/nio 


d5i exp 


r sj ] 

1 (71 

r ei 

1 -1 


2(7j 

\27V 

2af 


2 


(A4) 


where ui = Here 1 refers both to the (arbitrary) point where the field is evaluated, as well as to the scale it was 

smoothed with, R\. 

Next we evaluate the asymptotic behaviour at large 5cr for the joint two-point probability 

1 ( 5 ict 2 + S2cj'i — 2 ^{ri 2 )Si 62 
' 2 ' 


P(l,2) = 


1 


2TT^yaja1 -C(ri 2 ) Ac 


ddi / d (52 exp 

Acr 


- ?^(»’i 2 ) 


(A5) 


paying attention to the prefactors to the exponential terms. In general ci 7 ^ ( 72 , but when di and 62 represent the same field 
smoothed with the same filter taken at two different points (the case that we mostly need in this paper) cri = CT 2 . Introducing 
uncorrelated variables x = <^ig 2 +^ 2 gi j y _ obtain 

Oct(^2+Oct(^1 ^ ’ 


P(l,2) = 


(^cr(cri + CJ' 2 ) 

-^^( 712 ) A 


dx exp 


1 5cr (<7l -|- 1 X 2 )^ 
4criCT2 0-1(72 -|-C(^ 12 )' 


, ScT [x(cri-t-o-2)-2o-i] 

0 - 10-2 


f (5cr[2cr2—3:(crj^-fg2)] 
<Tlcr 2 


dy exp 


0-1(72 


4 (71(72 — C(^12) ‘ 


• (A 6 ) 


This integral is of Laplace type I = dx e f{x), a = which for large a asymtotically accumulates at the lower 
integration boundary over the interval x € [1,1 + Ax], Ax = ■ Asymtotic expansion is straightforward if one 

can expand /(x) in the Taylor series near x = 1 


f{x) = 


^ 5ci[x{ai+a2) — '2cri] 

^ 1^2 


(5cr[2g2-^(^l+^2)] 

o-lcr 2 




2(^_l)Moi_^exp 


(71(72 


1 e ((71 -(72) = 


4 (71(72 (71(72 — C(^12) 


in which case we get 


[ 0 - 10-2 + C(^ 12 )] 


3/2 


TT (5cr ((71 + (72)^((7i(72 — ^(ri2))3/2 


exp 


_l_c2 (^1 + ~ 2 ^( 712 ) 

2 “ (7=(7i-e=(ri2) 


(A7) 


(A 8 ) 


Two things are notable. First is the prefactor ((7i(72)/5cr- Second, we find that for small correlations the effect in the 
exponent where small ^l(J \02 is multiplied by (5cr/(7i(72 dominates the correction from the prefactor. Thus, as a leading order 
approximation, we can account for small correlations by factoring out the exponential correlation term from the original 
expression, with the values of the field replaced by the threshold values. 

In reality, the asymptotics in eq. do not give an accurate approximation if correlations are strong ^( 712 ) —> (7i(72, 

especially since our threshold parameter dcr/cr may not be very large. This is definitely the case for a distribution of identical 
objects at short distances, since then ^(r ^ 0) = ct=. More accurately, the Taylor expansion of /(x) in eq. EH is not 
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suitable when the width of the relevant integration range Aj/ = ScrAx{ai + o-i) j{axap) exceeds the width of the Gaussian 
\J\o\ai — C(ri 2 )] /(cri(T 2 ). In this case, however, the integration over y can be extended to ±oo. With subsequent asymptotic 
analysis of the integral over x this gives 


P(l,2) 


1 \Jo\ai + ^(ri2) 

A <5 cr (Tl + 0-2 


exp 


1 '^cr (l^l + 1^2)^ 

4 (71(72 (71(72 + ^(ri2) 


(71(72 + C(^12) 


> 


y/ (71(72 — ^(’" 12 ) vAA 


> 1 (A9) 


The general eqs. lEU-lElIl are much simpler in the case when variances are identical erf ~ cr^ ~ C(0)- Defining the cross 
correlation coefficient c(r 12 ) = C(’"i 2 )/C( 0 ) and specifying accurately the range of validity of eq. (EUl gives 


pa, 2) 
pa,2) 
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At 


V ^ A [c(ri 2 )] exp 


V ^ B [c(r 12 )] exp 


1 + c(ri2) 

,.2 


V ^ max 1, 


1 + c(r 12 ) I 

A/l-c(ri2) j ’ 


1 + c(ri2) 
l + c(ri 2 )J ’ - c(ri 2 ) 


> ^7 > 1 


(AlO) 

(All) 


where smooth functions A{x) = 1 and B(x) = > 1- As expected, p(l,2) p(l)p(2) and p(l,2) 

p(l). It is important to note that the probability is additionally enhanced by Z 2 = (5cr/(7 when correlations are strong. 
Finally, we combine equations ifvinii - ifrrni into the uniform approximation 


13(1,2) ^ C[c(ri2),i2]exp 


» 1, 


1 + c(ri2) 

with the help of an interpolating function, C(a;, v)^ such that C(0, v') = 1, C(l, v) = The choice 


(A12) 


C(a;, 12 ) 


v\px\/ (1 + xY 

(vy/ir — 1 ) VI — a; + (1 + a;) 


(A13) 


reflects both the details of the functions A{x) and B{x) and of the transition between equations ifTinii - ifrrni . 

In the weak correlation regime, the formula ifvr^ coincides with the classic result of Kaiser (1984). At the same time 
in the strong correlation regime, the result ifm shows that the correlation between regions of high density is additionally 
enhanced by r/^v factor. Although our result is rigorous for the points of high excursions of the field at all separations r, 
the interpretation of the last regime in terms of peak, or object, correlation, is questionable at r < + i ?2 when the two 

high-density points likely belong to the same peak. 
































